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INTRODUCTION 

As  part  of  a  study  of  the  distribution  of 
north  Pacific  albacore,  the  staff  of  the  Bureau 
of  Commercial  Fisheries  Biological  Labora- 
tory, Honolulu,  published  monthly  sea  surface 
temperature  charts  of  the  eastern  Pacific 
Ocean  from  1957  through  1959.  At  the  end  of 
this  phase  of  the  north  Pacific  albacore  investi- 
gations, the  task  of  preparing  and  distributing 
the  monthly  charts  was  assumed  by  the  staff 
of  the  Tuna  Forecasting  Program,  Bureau  of 
Commercial  Fisheries  Biological  Laboratory, 
San  Diego,  Calif. 

The  purpose  of  the  sea  temperature  charts 
now  is  to  provide  current  sea  temperature 
information  to  the  tuna  fishing  industry,  marine 
scientists,  and  meteorologists.  That  the  charts 
are  providing  useful  data  to  these  groups  is 
evident  from  the  work  of  Broadhead  and 
Barrett  (1964),  Dickson  (1964),  Flittner  (1963), 
Magill  (1963),  Namias  (1959,  1963),  Quast 
(1964),  and  others. 

Charts  prepared  at  the  BCF  Biological 
Laboratory  at  San  Diego  are  circulated  in  two 
series.  Monthly  sea  temperature  charts  are 
published  in  the  Bureau  of  Commercial  Fish- 
eries California  Fishery  Market  News  Monthly 
Summary,  Part  II  -  Fishing  Information.  A 
15-day  temperature  chart  for  use  by  albacore 
fishermen  is  published  from  April  through 
October  as  a  supplement  to  the  above  publica- 
tion for  a  limited  area  of  the  temperate  eastern 
Pacific.  Monthly  summaries  and  charts  of 
meteorological  variables  and  heat  budget  com- 
putations are  circulated  to  a  limited  mailing 
list  under  the  title  "Summary  of  Marine  Synop- 
tic Weather  Reports." 

The  purposes  of  this  paper  are  to  describe 
an  automatic  data  processing  program  (ADP) 
for  summarization  of  synoptic  marine  weather 
reports,  to  describe  the  sea  temperature 
reporting  program,  and  to  provide  information 
on  procedures  for  calculating  the  energy 
exchange  at  the  air-sea  interface  for  the 
Pacific  Ocean. 


Note. — James  H.  Johnson,  Oceanographer;  Glenn  A. 
Flittner,  Fishery  Biologist  (Research);  and  Marvin  W. 
Cline,  Mathematics  Technician,  Bureau  of  Commer- 
cial Fisheries  Biological  Laboratory,  San  Diego, 
Calif. 


SOURCES  OF  DATA 

The  first  International  Maritime  Conference, 
called  in  Brussels,  Belgium,  in  1853,  conceived 
the  idea  that  cooperating  governments  should 
encourage  routine  weather  observations  on 
ships,  and  that  nations  should  prepare  and 
publish  charts  of  ocean  currents,  prevailing 
winds,  average  sea  and  air  temperatures,  and 
storm  tracks  based  on  these  and  other  data. 
Since  inception  of  the  program,  weather  report- 
ing has  played  an  increasingly  vital  role  in 
the  safe  and  efficient  operation  of  ships 
at  sea. 

Data  used  in  preparing  the  charts  and 
summaries  come  from  a  variety  of  sources. 
The  main  source  is  the  synoptic  marine  radio 
weather  report  from  ships  at  sea.  Additional 
data  are  provided  by  organizations  that  mail 
coastal  and  ocean  temperatures  to  the  San 
Diego  laboratory  for  inclusion  in  the  charts. 
Early  in  the  sea  temperature  reporting  pro- 
gram, it  became  apparent  that  use  of  ADP 
would  be  helpful  in  preliminary  data  screening 
for  radio  transmission  errors,  for  coding 
errors,  and  for  compiling  an  average  of  about 
10,000-15,000  synoptic  marine  weather  obser- 
vations each  month.  After  the  feasibility  of 
summarizing  sea  temperature  data  by  ADP 
methods  was  established,  it  was  expedient  to 
broaden  the  program  to  compile  other  meteoro- 
logical variables  for  use  by  scientists  in 
marine  climatology  and  other  research  studies. 
The  elements  for  computation  of  the  energy 
exchange  at  the  air-se'a  interface  were  avail- 
able from  these  summaries  and  monthly  aver- 
ages of  this  exchange  by  5-degree  quadrangles 
are  now  computed  routinely. 

The  basic  data  utilized  by  the  system 
described  below  is  the  "International  Surface 
Report  from  Ship  in  Full  Form,  FM  21.A," 
approved  by  the  World  Meteorological  Organi- 
zation (WMO).  Cooperating  American  and 
foreign-flag  vessels  make,  record,  and  trans- 
mit the  standard  marine  weather  observations 
according  to  established  procedures  set  up  by 
WMO.  Observations  taken  at  0000,  0600,  1200, 
and  1800  Greenwich  Mean  Time  daily  are 
transmitted  to  designated  commercial  and 
Government  radio  stations  around  the  world 
(see  U.S.  Navy  Hydrographic  Office  Publication 
No.  206,  Radio  Weather  Aids,  1958). 


Depending  on  geographical  location,  ships 
at  sea  report  marine  weather  to  designated 
receiving  stations.  North  Pacific  Ocean  data 
for  the  region  west  of  the  180th  meridian  are 
received  by  "Observer  Tokyo."  Data  for  the 
region  east  of  the  180th  meridian  are  received 
by  "Observer  San  Francisco"  and  are  relayed 
to  the  U.S.  Weather  Bureau  Forecast  Center 
at  San  Francisco  International  Airport,  San 
Bruno,  Calif.  This  center  transmits  both 
Observer  San  Francisco  and  Tokyo  reports  to 
other  agencies  via  teletype  circuits.  The  BCF 
Biological  Laboratory,  San  Diego,  maintains 
a  drop  on  Service  "O",  Circuit  8274  to  receive 
ship  surface  weather  information  24  hours 
daily. 

The  ship  weather  report  code  adheres  to 
the  following  symbolic  5-digit  word  group 
format: 


Table  1. — Elements  of  the  ship's  surface 
synoptic  code,  FM  21. A1 


Symbols 


Definitions  of  symbols 


YQLaLaLa 
PPPTT 


OT    T    T,T, 
s     s    d    d 


L   L   L   GG 

o    o    o 


Nddff       WwwW 


N,CT  hC.  ,C__ 
h    L       M    H 

Id    d    P    H 

WWWW 


Dsvsapp 

ICE       c,KD.re 
2       i 


Additional  groups  are  coded  when  required; 
supplementary  information  following  the  ICE 
word  group  is  usually  transmitted  in  plain 
language.  Definitions  of  symbols  are  given  in 
table   1. 

Since  a  large  percentage  of  the  incoming 
marine  weather  observations  emanates  from 
vessels  having  observers  with  limited  meteor- 
ological experience,  the  question  arises  as 
to  the  validity  of  certain  portions  of  the  stand- 
ard weather  message.  It  is  recognized  that 
certain  observational  errors  maybe  introduced 
by  equipment  in  poor  condition  or  through 
faulty  operation.  The  problem  of  variability  of 
ship's  injection  system  thermometer  readings, 
as  an  example,  has  been  reviewed  by  Saur 
(1963).  From  his  studies  aboard  U.S.  Navy 
radar  picket  vessels,  Military  Sea  Transport 
Service  ships,  and  oil  tankers,  Saur  found  that 
the  average  injection  temperature  bias  from 
that  of  a  bucket  temperature  taken  at  the 
surface  was  about  +  1.2°  F.  In  addition  to 
instrument  error,  one  complicating  factor  is 
that  sea  temperatures  are  not  reported  directly 
but,  rather  must  be  derived  from  air  tempera- 
tures and  air-sea  temperature  differences.  In 
the  near  future  the  Commission  on  Maritime 
Meteorology  of  the  WMO  plans  to  discuss  the 
possibility  of  changing  the  marine  weather 
coding  instructions  so  that  sea  temperature  will 
be  reported  directly.  Another  complicating 
factor  is  the  problem  of  communication 
errors.  A  study  by  Gibson  (1962)  for  the 
Atlantic  Ocean  showed  that  19  percent  of 
the  sea  temperatures  received  on  1  day 
were  incorrect;  these  errors  were  due  to 
improper  coding  procedures,  to  communica- 
tion transmission  errors,  and  to  faulty  in- 
struments. 


Y 
Q 
LaLaLa 
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N 
dd 

ff 

W 

WW 

W 
PPP 

TT 
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CL 
h 
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PP 


Vs 

Vd 


1 

<Vw 

pw 

«w 

ICE 

c2 
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Di 

r 
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Day  of  week 

Octant  of  globe 

Latitude  in  degrees  and  tenths 

Longitude  in  degrees  and  tenths 
Greenwich  civil  time  of  observa- 
tions,  to  nearest  hour 

Fraction  of  celestial  dome  covered 

by  clouds 
Direction  from  which  the  wind  is 

blowing  (tens  of  degrees) 
Wind  speed  in  knots 

Visibility  (in  code) 
Present  weather  (in  code) 
Past  weather  (in  code) 

Sea-level  pressure   (tens,  units, 

tenths)   millibars   (mb). 
Temperature  of  the  air  (whole 

degrees  Centigrade) 

Fraction  of  celestial  dome  covered 
CL  or  ^U  cloud 

Clouds  of  types  stratocumulus ,   stra- 
tus, cumulus,    cumulonimbus  (in  code) 

Height  of  base  of  lowest  C^  or  Cjj 
cloud  above  sea 

Clouds  of  types  altocumulus,   alto- 
stratus,  nimbostratus  (in  code) 

Clouds  of  types  cirrus,   cirrostratus, 
cirrocumulus  (in  code) 

Ship's  course   (in  code) 
Speed  of  ship  in  knots  (in  code) 
Characteristic  of  barometric  tend- 
ency (in  code) 
Amount  of  barometric  change  (units 
and  tenths)  millibars  (mb.) 

Group  designator 

Difference  between  air  and  sea  tem- 
perature (half  degrees  Centigrade) 

Temperature  of  dew  point  (whole 
degrees  Centigrade) 

Group  designator 

Direction  from  which  waves  are 

coming  (tens  of  degrees) 
Period  of  waves   (in  code) 
I* an  height  of  waves  (in  code) 

Group  designator  for  ice  group 
Description  of  kind  of  ice  (in  code) 
Effect  of  ice  on  navigation  (.'.a  code) 
Bearing  of  ice  limit  (in  code) 
Distance  of  ice  from  ship  (in  code) 
Orientation  of  ice  field   (in  code) 


1  Adapted  from  U.S.   Weather  Bureau  Manual  of 
Marine  Meteorological  Observations,   Circular  M 


[11th  ed.),   January  1963,  p.   8. 


Errors  introduced  into  the  weather  message 
for  various  reasons  are  often  difficult  to 
note  unless  they  are  of  serious  nature.  Data 
screening  at  the  BCF  Biological  Laboratory 
in  San  Diego  is  designed  primarily  to  detect 
the  following  errors: 

1.  Departures  from  major  5-digit  group- 
ings: e.g.,  teletype  garbles,  run-on  groups 
of  5  or  more  digits  not  separated  by  a  space, 
failure  to  shift  from  lower  to  upper  case,  and 
vice  versa. 

2.  Major  coding  errors  inbarometric  pres- 
sure readings:  all  coded  messages  are  inter- 
preted for  values  between  950.1  and  1049.9 
mb.;  coded  values  greater  than  500  are  assigned 
to  the  range  950.1-999.9  mb.,  and  values  less 
than  500  are  assigned  to  the  range  1000.0- 
1049.9  mb.  Observed  barometric  pressures 
will  be  refined  further  in  the  future  when 
monthly  barometric  averages  and  associated 
statistics  become  available  for  screening 
purposes. 

3.  Major  errors  in  sea  temperatures:  sea 
temperature  observations  are  screened  by 
comparison  with  a  table  of  values  interpolated 
from  monthly  means  for  the  north  Pacific 
Ocean  found  in  the  "Atlas  of  Climatic  Charts 
of  the  Oceans,"  U.S.  Weather  Bureau.  Tem- 
peratures outside  a  selected  range  are  rejected 
and  placed  in  a  summary  listing  for  inspection 
by  the  analyst. 


DISTRIBUTION  OF  OBSERVATIONS 

Because  the  data  used  are  compiled  pri- 
marily from  merchant  marine  weather  reports, 
the  distribution  of  observations  over  the  ocean 
is  irregular.  Observations  are  most  numerous 
in  the  major  shipping  lanes  from  San  Fran- 
cisco to  Hawaii,  San  Francisco  to  Japan, 
Panama  to  San  Diego,  Los  Angeles  and  San 
Francisco,  and  from  Seattle  to  Japan.  Shipping 
lanes  of  secondary  importance  are  from  San 
Diego  to  Hawaii,  Panama  to  Hawaii,  and 
Panama  south  to  the  ports  off  the  west  coast 
of  South  America.  Vast  regions  of  the  equa- 
torial zone  and  south  Pacific  Ocean  have 
much  less  merchant  vessel  traffic  than  the 
north  Pacific  and,  consequently,  fewer  reports 
are  received  from  these  areas.  Sample  dis- 
tributions of  the  number  of  usable  sea  tem- 
perature observations  are  given  in  figures 
1  and  2. 

In  recent  years,  progress  has  been  made 
in  obtaining  coverage  from  some  regions  that 
heretofore  have  yielded  few  reports.  One  quite 
important  region  that  has  been  little  covered 
is  the  area  off  the  west  coast  of  South  America. 
The  importance  of  this  area  canbe  appreciated 
>y     examining     the     world    fishing    statistics. 


In  1961,  FAO  reported  that  Peru  was  surpassed 
only  by  Japan  in  world  production  of  fish 
and  fish  byproducts.  Recently,  WMO  has  asked 
that  cooperating  ships  steaming  in  the  area 
south  of  the  equator  and  from  the  west  coast 
of  South  America  to  long.  140°  W.  send  reports 
to  "Observer  San  Francisco."  Furthermore, 
the  U.S.  Weather  Bureau  has  recently  estab- 
lished a  Port  Visitation  Office  in  Panama  and 
is  encouraging  vessels  passing  through  the 
canal  into  the  Pacific  to  send  reports  to 
designated  receiving  agencies. 


AUTOMATIC  DATA  PROCESSING  OF 
WEATHER  REPORTS 

The  Program 

Incoming  data  are  received  from  the  U.S. 
Weather  Bureau,  San  Francisco  International 
Airport,  San  Bruno,  Calif.,  via  the  landline 
teletype  circuit,  Service  "O",  Circuit  8274, 
mentioned  previously.  Each  transmission  is 
converted  to  a  5-channel  paper  tape  by  a 
teletype  reperforator  receiver.  Concurrent 
with  the  reperforator  operation,  the  trans- 
mission is  duplicated  on  a  page  printer 
machine  for  message  verification  purposes 
and  as  a  backup  when  the  reperforator  mal- 
functions. 

The  perforated  tape  is  fed  into  an  IBM 
047  tape-to-card  key  punch  machine  to  convert 
each  ship's  weather  observation  to  standard 
data  card  format.  These  data  cards  are  then 
accumulated  each  month  for  processing.  The 
data  are  presorted  by  octant  and  longitude 
into  10-degree  wide  strips  to  facilitate  proc- 
essing by  the  machine  and  to  provide  a 
logical  output  format.  The  information  on 
each  of  the  sorted  cards  is  then  loaded  onto 
magnetic  tape  for  high-speed  input  to  the 
computer.  The  basic  steps  are  schematized  in 
figure  3. 

The  program  in  current  use  is  written  in 
FORTRAN-62  source  language  for  the  Control 
Data  Corporation  1604  computer  located  at 
the  University  of  California,  San  Diego.  A 
machine  language  deck  has  been  prepared  from 
the  source  language  statements  for  routine 
processing  of  synoptic  data.  The  computer 
program  processes  the  cards  and  prepares 
output  summaries  for  each  10-degree  strip  of 
longitude  separately  until  all  accumulated 
cards  are  processed.  Data  cards  are  sum- 
marized at  the  rate  of  5  per  second  and  the 
output  tape  is  produced  at  the  rate  of  40 
lines  per  second.  The  output  tape  is  subse- 
quently listed  on  multicopy  paper  for  tabulation 
assembly.  The  flow  charts  in  figures  4-6 
depict  the  major  data  processing  functions 
performed;  numerical  designators  are  keyed 
to  processing  statements  within  the  program 
listing  in  Appendix  A. 
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Figure  1. — Distribution  of  screened  sea  temperature  observations  derived  from  the  synoptic  marine  weather  observations, 

August  1963. 
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Figure  2. — Distribution  of  screened  sea  temperature  observations  derived  from  the  synoptic  marine  weather  observations, 

February  1964. 
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Figure    3.— Schematic   diagram   of   the   automatic   weather  data   processing   system   in   use  by  the  Bureau  of  Commercial 

Fisheries  Biological  Laboratory,  San  Diego,  Calif. 
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Figure  4.— Schematic  flow  chart  depicting  the  major  data  processing  functions  performed  by  the  Control  Data  Corporation 
1604  computer  at  the  University  of  California,  San  Diego.  See  Appendix  A  for  key  to  numerical  designations. 
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Figure  5.— Schematic  flow  chart  depicting  the  major  data  processing  functions  performed  by  the  Control  Data 
Corporation  1604  computer  at  the  University  of  California,  San  Diego.  See  Appendix  A  for  key  to  numerical 
designations.  (Continued) 
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1604  computer  at  the  University  of  California,  San  Diego.  See  Appendix  A  for  key  to  numerical  designations.  (Continued) 


Monthly  Data  Readout 

It  is  becoming  increasingly  clear  to  marine 
scientists  that  much  interaction  occurs  between 
the  ocean  and  the  atmosphere  (Eber  and  Sette, 
1959;  Hanzawa,  1962;  Namias,  1963;  Uda,  1960, 
1962).  To  provide  data  for  research  into  air-sea 
interaction  studies,  not  only  sea  surface  tem- 
peratures but  also  other  variables  contained  in 
the  marine  synoptic  weather  report  are  sum- 
marized. The  basic  elements  of  the  monthly 
data  compilations  are  listed  in  four  major 
tables,  as  follows1: 

Table  1.001-  1.069--Sea  surface  temperatures 
(°  F.)  with  number  of  observations  by  1-degree 
quadrangles  for  the  month  of  [entry  changes 
each  month]. 

Table  2.001-2.184--Barometric  pressures 
(mb.),  X-  and  Y-wind  vectors  (knots),  and 
average  wind  speed  (knots)  with  number  of 
observations  by  1-degree  quadrangles  for  the 
month  of  [same  as  above]. 


Table  3.001-3.046--Sea  surface  temperatures 
(  F.  and  C.)  and  number  of  observations  by 
2-degree  quadrangles  for  the  month  of  [same 
as  above]. 


xFor  convenience  in  output  assembly,  pages  within  each 
table  are  numbered  consecutively  from  .001. 


Table  4.001-4.069--Sea  surface  temperatures 
(  C),  air  temperatures  (°  C.),  dew  point 
temperatures  (°  C),  barometric  pressures 
(mb.),  X-  and  Y-wind  vectors  (knots),  average 
wind  speeds  (knots),  cloud  cover  (tenths  of 
sky  covered),  and  heat  budget  (calories  per 
square  centimeter  per  day),  with  number  of 
observations  by  5-degree  quadrangles  for  the 
month  of  [same  as  above]. 

Sample  page  copy  of  each  table  listed  above  is 
included  in  Appendix  C.  Each  table  conforms 
to  the  following  conventions: 

1.  Latitude  and  longitude  range  headings -- 
Identify  a  10-degree  wide  strip  of  longitude 
and  give  its  length  in  latitude.  Each  successive 
page  with  identical  range  headings  belongs  in 
the  same  10-degree  strip. 

2.  Latitude  and  longitude  indices--Give 
map  location   of  each  piece  of  data.  Answer  is 


assumed  to  be  for  the  geometric  center  of 
each  quadrangle.  Location  given  on  output  is 
actually  the  corner  of  the  quadrangle  and  fol- 
lows this  convention: 

Octants  0,  1  -  lower  right  corner 

Octants  2,  3  -  lower  left  corner 

Octants  5,  6  -  upper  right  corner 

Octants  7,  8  -  upper  left  corner 

(Where  octant  designations  conform  to  Hydro- 
graphic  Office  Publication  No.  206  code  form 
FM  21.A.  "Surface  Report  From  Ship  in  Full 
Form.") 

3.  Wind   vectors-- Wind   vectors    follow  this 
convention: 

(a)  X--vectors  run  due  east  and  west 
East--positive 
West--negative 

(b)  Y--vectors  run  due  north  and  south 
North- -positive 

South-  -negative 

4.  Average  wind  speed--The  simple  average 
of    all   observations    irrespective    of  direction 


Y2!. 


(Note  -  does  not  necessarily  equal 

5.  Heat  budget  computations --The  net 
energy  exchange  at  the  air-sea  interface.  The 
printout  has  the  following  sign  convention: 

(a)  Energy  into  sea--positive 

(b)  Energy  out  of  sea--negative 

The     computer    printout    labels    of    various 
components  of  the  equation  are: 

Q(T)  -  total  (net)  energy  exchange  across 
air-sea  interface) 

Q(  I )  -  incoming  solar  radiation  corrected 
for  cloud  cover 

Q(R)    -    radiation  reflected  by  sea  surface 

Q(B)    -    back  radiation 

Q(E)  -  heat  loss  or  gain  through  evapora- 
tion or  condensation 


Q(H)    -    sensible   heat   conduction   between 
sea  and  air 

Energy  exchange  equations  will  be  discussed 
in  subsequent  sections  of  this  paper. 

The  following  notes  are  given  to  clarify 
some  aspects  of  the  data: 

1.  A  blank  in  heat  budget  output  [Q(T),  Q(I), 
Q(R),  Q(B),  Q(E),  Q(H)]  is  equal  to  0. 

2.  When  any  heat  budget  value  [Q(T),  Q(I), 
Q(R),  Q(B),  Q(E),  Q(H)]  is  between  -1  and  0, 
a  minus  sign  only  appears  as  entry. 

3.  Any  entry  or  a  blank  followed  by  (0) 
means  no  data  for  that  quadrangle. 


CHARTS  DERIVED  FROM 
MONTHLY  READOUT 

Sea  Surface  Temperature  Charts 

Examples  of  monthly  sea  surface  tempera- 
ture charts  obtained  by  using  the  data 
summaries  and  reported  in  the  Bureau  of 
Commercial  Fisheries  California  Fishery 
Market  News  Monthly  Summary,  Part  II  - 
Fishing  Information  are  shown  in  figures 
7-12.  Figure  7  depicts  the  monthly  sea  surface 
isotherm  configuration  derived  from  ships' 
weather  reports,  average  sea  temperatures 
at  Weather  Ship  Stations  November  and  Papa, 
and  average  temperatures  of  coastal  stations. 
Because  coastal  station  data  are  mailed  and 
must  reach  the  BCF  San  Diego  Laboratory  by 
the  first  day  after  the  close  of  a  month  for  proc- 
essing, the  coastal  station  averages,  except  for 
that  reported  at  the  pier  of  Scripps  Institution 
of  Oceanography,  usually  represent  only  about 
the  first  25  days  of  each  month.  Coastal  sta- 
tions included  in  figure  7  are  as  follows: 


Station 


Location 


Latitude 


Longitude 


Langara  Island 54°  15 

Cape  St.  James 51°56 

Quatsino  Island 50°27 

Amphitrite  Point 48°55 

Umatilla  Lightship 48°10 

Columbia  River  Lightship 46°  11 

Blunt' s  Reef  Lightship 40°26 

Russian  Gulch  State  Park 39°18 

Sonoma  Coast  State  Park 38°23 

Farallon  Island  Lightship 37°46 

Hopkins  Marine  Station 36°38 

Gaviota  State  Park - 34°28 

Santa  Barbara  -- 34°24 

Doheny  State  Park - --  33°28 

Scripps    Institution  of  Oceanography 32°52 

Ens.enada - - 31°52 

Guadalupe  Island -  28°41 

Cedros  Island - 28°15 


N. 
N. 

N. 
N. 

N. 
N. 

N. 

N. 
N. 
N. 
N. 
N. 
N. 
N. 
N. 
N. 
N. 
N. 


133°03 
131°01 
128°02 
125°32 
124°50 
124°  11 
124°30 
123°48 
123°05 
122°42 
121°55 
120°14 
119°42 
117°41 
117°15 
116°38 
118°17 
115°10 


W. 
W. 

w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
w. 
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Figure    11. —Deviation   of   sea   surface  temperatures  (°  F.)  from  February  1963  -  February  1964.  Shaded 

areas  colder  in  1964. 
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Figure    12. — Etevlatlon    of    sea    surface    temperatures   (°  F.)   from    12-year  mean  (1947-58)  by  2-degree 
quadrangles,  February  1964.  Shaded  areas  colder  In  1964. 
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Monthly  sea  surface  temperature  deviations 
for  February  1964,  from  the  long-term  mean 
charts  given  in  H.O.  225,  World  Atlas  of  Sea 
Surface  Temperatures,  2d  ed.,  1944,  appear 
in  figure  8.  The  exact  base  period  used  for  con- 
struction of  the  mean  charts  is  not  known.2 

The  deviation  of  sea  surface  temperature, 
February  1964,  from  the  corresponding  month 
of  the  previous  year  is  given  in  figure  9. 

Whereas  three  of  the  charts  (figs.  7-9) 
encompass  the  temperate  north  Pacific  Ocean 
east  of  the  180th  meridian,  the  remainder 
(figs.  10-12)  give  data  for  the  eastern  tropical 
Pacific  Ocean.  The  first  (fig.  10)  of  these 
gives  actual  average  temperatures  by  1 -degree 
quadrangles;  underscored  values  represent 
averages  from  two  or  more  observations. 
The  next  chart  (fig.  11)  gives  the  deviation 
from  the  previous  year  by  1 -degree  quad- 
rangles, and  the  last  (fig.  12)  shows  the 
deviation  by  2-degree  blocks  from  the  12-year 
(1947-58)  mean  prepared  by  Renner  (1963). 
When  the  10-year  (1950-59)  mean  is  completed 
by  the  BCF  Biological  Laboratory  at  Stanford, 
this  chart,  too,  will  be  based  on  the  same 
period.  Because  of  the  scarcity  of  observations 
from  the  tropical  regions,  actual  values  are 
presented  instead  of  contoured  values. 

From  April  through  October,  as  a  supple- 
ment to  the  Bureau  of  Commercial  Fisheries 
California  Fishery  Market  News  Monthly  Sum- 
mary, Part  II  -  Fishing  Information,  a  series 
of  charts  covering  15-day  intervals  is  pre- 
pared and  distributed  to  albacore  fishermen 
and  to  others  interested  in  receiving  it  (fig.  13). 
This  chart  covers  the  region  from  the  coast 
of  North  America  to  long.  136°30'  W.andfrom 
lat.  25°  N.  to  52°  N.  It  is  designed  to  span  the 
area  and  season  when  the  commercial  albacore 
fishery  takes  place  off  the  Pacific  coast. 

Meteorological  Charts 

Charts  are  prepared  each  month  from  the 
readout  for  averages  of  cloud  cover  (figs.  14 
and  15),  wind  speed  (figs.  16  and  17),  air 
temperature  (figs.  18  and  19),  dew  point 
temperature  (figs.  20  and  21),  sea  temperature 
(figs.  22  and  23),  and  barometric  pressure 
(figs.  24  and  25).  Values  are  summarized  by 
5-degree  quadrangles  of  latitude  and  longitude. 


2 The  Introduction  of  H.O.  225  states  as  follows:  "This 
world  atlas  of  sea  surface  temperatures  was  prepared  by 
the  Scripps  Institution  of  Oceanography,  La  Jolla,  Cali- 
fornia, from  all  sea  surface  temperature  information 
collected  by  the  U.S.  Hydrographic  Office,  together  with 
all  available  data  from  other  sources,  with  special 
reference  to  regions  of  few  observations.  It  is  believed 
that  this  compilation  represents  the  most  accurate  pattern 
of  sea  surface  temperatures  of  the  world  yet  prepared." 
The  staff  of  the  Bureau  of  Commercial  Fisheries  Biological 
Laboratory,  Stanford,  Calif.,  is  now  preparing  a  10-year 
mean  for  the  period  1950-59.  When  the  newmean  becomes 
available,  the  anomaly  chart  will  be  based  on  it. 


Values  are  not  plotted  for  5-degree  quad- 
rangles having  fewer  than  five  observations. 
Two  charts  of  each  of  the  variables  at  opposite 
climatological  seasons  show  the  magnitude 
of  seasonal  changes  occurring  over  the  ocean. 
Charts  for  August  1963  cover  the  Pacific 
only  from  the  American  coasts  to  long.  180° 
because  it  was  not  until  January  1964  that  data 
were  programmed  to  cover  the  entire  Pacific 
Ocean.  These  data  are  distributed  monthly  on 
a  limited  mailing  list. 

Energy  Exchange  Computations  and  Charts 

The  marine  synoptic  weather  report  contains 
basic  data  for  computation  of  the  energy  ex- 
change occurring  at  the  air-sea  interface. 
The  methods  described  below  represent  one 
of  the  first  attempts  to  monitor  the  gross 
energy  flux  at  the  air-sea  interface  over  a 
broad  expanse  of  the  ocean  on  a  monthly  basis. 

For  a  given  area  and  time  period,  the 
equation  for  the  energy  exchange  at  the  air- 
sea  interface  is  QT  =   Qj  -  QR  -  Qfi  -  Q£  -  Q 

(See  page  10  for  definitions  of  elements  of  this 
equation.)  The  gainor  loss  of  heat,  added  to  that 
resulting  from  horizontal  and  vertical  advec- 
tion,  represents  the  total  heat  energy  available 
to  the  upper  mixed  layer  of  the  ocean.  Energy 
exchange  calculations  presented  here  do  not 
take  into  account  changes  in  heat  brought 
about  by  advection. 

Energy  exchange  equations  have  been  re- 
viewed in  detail  by  Sverdrup,  Johnson,  and 
Fleming  (1942),  Jacobs  (1951),  Roden(1959), 
Laevastu  (I960),  and  others.  Opinion  differs 
considerably  as  to  the  accuracy  of  the  various 
equations  describing  the  total  energy  exchange 
at  the  air-sea  interface.  It  is  not  the  intent 
here  to  review  the  accuracy  of  the  many 
different  equations  that  have  been  proposed. 
Controlled  experiments  similar  to  those  car- 
ried out  on  Lake  Hefner  by  Anderson  (1954)  and 
Marciano  and  Harbeck  (1954)  are  needed  to  de- 
termine which  of  the  equations  most  closely 
describes  what  actually  takes  place  in  the  ocean. 

Qt,  incoming  solar  radiation  corrected  for 
cloud  cover  (in  cal./cm. /day),  is  determined 
from  the  following  equation  proposed  by 
Berliand  (1960): 

Q:  =  (Berliand  table)3   x   (1  -  aC  -  bC2) 

where  C_  -  cloudiness  in  tenths 

Jb   =    0.38 
and         a  =    a  function  of  latitude  as  follows: 


Latitude: 


0 


10         15         20         25 


.a  0.38     0.40     0.40     0.39     0.37     0.35 

Latitude:         30        35        40        45        50        55 
a  0.36     0.38     0.38     0.38     0.40     0.41 


See  Appendix  table  B-l  for  Berliand's  table  which  lists 
values  of  incoming  solar  radiation  with  a  clear  sky. 
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Figure  13. — Sea  surface  temperature  chart  for  15-day  period,  August  1-15,  1963. 
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Figure  14. — Cloud  cover  (hundredths),  August  1963. 
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20 


8 

* 

1   f)  1    I 

J-H-Lf 

/13/191 

Jifi3\i4^ 

9\l? 

M^skijii 

12  I10I9  \  8  1 

b\V 

te                                            #>T 

'is    16  he 

17 

IS 

2  \lo\ll  \»\ 

bVI*                       V\ 

/|4  /  14  (1  3 

IS 

14 

12  110  \  9  \  9 

1  9  \  II  \8                                                                |  \    \j 

/14J  12     16 

13 

14 

II    1  8   \  9  \  8 

\  "All-w5^                          «<*r^S\\  \ 

/  II    /  12   |l2 

14 

II 

II   lio  1  II  \io  \  11 

ViViA^Sk            A   \  V5^V\ 

/  II    J  12    |  II 

II 

13 

14    |  14  |  12  \  16   \  10  \  3 

\  l2Al2-WwS;        A\\     yV^Mi  Jy 

/         J  12   |  14 

M 

13  j  '-5. 

15 

15J1S  |l5  \  14^10  \>0  \^_\^-\-"^S     'N^^fA,     W\ 

12 

Il 

1 

13 



14 

14 

IS  1 14  113  1  12  \ 

10  \ 

^Mi43j^s*V\\j\ 

13 

16 

13 

13 

8 

6 

13  111  1 

s  \7  \e  y 

a  \»\ 

4 

7 

10 

12 

8 

II 

7 

IS 

II 

12 

13  j  11 

9 

7 

9    I  10 

10  j  7 

e^sy 

V^d 

Srr 

140" 

150* 

160* 

17C 

180" 

17 

r 

Kcr  * 

ISO- 

13  mq. 

wr 

120- 

IKT 

"ocril 

II  i,J2 

10  956 

0 

/(A 

«•*     6-» 

y**-^ 

• 

X 

14 

13 

14 

8    le 

a/- 

*        I* 

&   ^ 

**       k^ 

l\. 

15 

14 

16 

/  (2    1 

-. 

— /     1     |     lB  IV       I 

/    /    /    /   tx    if 

-vTl   %   \     \     \ 

\   V-  I 

1     ^1 

1      1      t — / — T~~i~~ Z?l-     / 

t~H-4JJ  f~hJ 

/  /  7 1  1i~~T~~i-~J  /  /  iii°' 

!~i~f~~/~~/~jj  /  ~r~y^cy  /17*'- 

/7>^2//XX//^<^^ 

Vv^-^-^sr 

s                                       s. 
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Figure  17. — Wind  speed  (knots),  February  1964. 
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Figure  18. —Air  temperature  (°  C),  August  1963. 


23 


Figure  19. — Air  temperature  (°  C.),  February  1964. 
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Figure  20. — Dew  point  temperature  (°  C),  August  1963. 


25 


Figure  21. — Dew  point  temperature  (°  C),  February  1964. 
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Figure  22. — Sea  surface  temperature  (°  C),  August  1963. 
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Figure  23. — Sea  surface  temperature  (°  C),  February  1964. 
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Figure  24. — Barometric  pressure  (mb.),  August  1963. 
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Figure  25.  —  Barometric  pressure  (mb.),  February  1964. 
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Of  the  incoming  radiation  corrected  for  the 
screening  effects  of  cloud  cover,  some  is 
reflected  at  the  sea  surface.  The  amount 
reflected  depends  on  the  latitude  and  time  of 
year.  The  computational  expression  is  as 
follows: 


4.70  (e 
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Qxx 


where  r  =  percentage  of  radiation  reflected 
presented  in  a  table  by  Budyko  (1956);  this 
table  is  listed  in  Appendix  B- 2.  The  percentage 
varies  from  about  6  percent  in  low  latitudes 
to  more  than  20  percent  in  high  latitudes  in 
winter. 

Back  radiation,  Qg  (in  cal./cm.  /day),  is 
the  difference  between  long-wave  radiation 
from  the  sea  surface  and  long-wave  radiation 
from  the  atmosphere.  The  following  semi- 
empirical  equation  proposed  by  Berliand  and 
Berliand  (1952),  and  used  by  Roden  (1959)  for 
studies  of  the  California  Current  system,  has 
been  incorporated  in  this  study: 


=  so8 


(0.39    -  0.050V~i)  (1   -  k  C2) 


+  4so6 
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where     _s    =  0.97  -  the  ratio  of  the  radiation  of 

the  sea  surface  to  a  black  body 
8  8  =  absolute  temperature  in  °  C; 

ai£  =    1.175   x    10-7  .  the  Stefan-Boltz- 
mann  constant; 
e  =  vapor  pressure  (mb.); 
k  =  constant  *; 
C  =  cloudiness  in  tenths. 

Off  the  California  coast,  Roden  (1959)  has 
found    that    QR    varies    from    about   80   to    140 

2       — 

cal./cm. /day  and  Seckel  (1962)  has  found  that, 

in  the  Hawaiian  Islands  region,  values  range 
from  115  to  150  cal./cm2/day.  At  Triple 
Island,  British  Columbia,  and  at  ocean  station 
"Papa,"  Tabata  (1958,  1961)  has  recorded 
average  values  of  about  100  cal./cmr/day 
with  little  seasonal  variation. 

The  incoming  radiation  corrected  for  cloud 
cover  minus  that  which  is  reflected  and  lost 
by  back  radiation  may  be  called  the  "effective" 
radiation.  Additional  energy  enters  or  leaves 
the  sea  surface  as  evaporation  (Qg)  and  sensi- 
ble heat  (Qu),  Jacobs  (1951)  discussed  evapora- 
tion and  conduction  of  sensible  heat  at  length 
and  presented  seasonal  charts  of  these  values 
for  the  north  Pacific  Ocean. 

Evaporation  depends  upon  (1)  the  velocity 
of  the  wind  and  of  the  vapor  pressure  differ- 
ence between  the  sea  surface  and  air  above  it, 
and  (2)  a  constant  coefficient  of  proportionality. 
There  is  little  agreement  as  to  the  value  of 
the  constant  that  should  be  used;  the  one 
selected  for  our  computations  is  given  by 
Tabata  (1958): 


where  e„  =  vapor   pressure   at   temperature   of 
sea  surface  (mb.); 
ea  =  vapor  pressure  of  air  (mb.); 
and     W     =  wind  speed  (m./sec. ). 

Appendix  Table  B-4  lists  values  of  satura- 
tion vapor  pressure  over  water  used  in  compu- 
tation of  evaporation. 

In  general,  regions  of  greatest  evaporation 
are  those  wherein  northerly  transport  of 
surface  water  is  greatest  and  which  are 
subjected  during  winter  to  frequent  invasions 
by  cold,  dry  air  masses  from  the  interiors 
of  continents.  In  the  Pacific,  this  region  is 
situated  off  Japan  where  cold  dry  air  of 
continental  Asiatic  origin  frequently  traverses 
the  northward-flowing,  warm  Kuroshio  Cur- 
rent. Another  region  of  high  evaporation  is 
in  the  trade  wind  zone  as  a  result  of  relatively 
strong  winds  and  dry,  descending  air  associated 
with  the  semipermanent  fields  of  high  pressure 
(Jacobs,  1951).  A  region  of  low  evaporation, 
in  contrast,  is  in  the  eastern  Pacific  over  the 
southward-flowing,  cool  California  Current. 

Bowen  (1926)  established  the  relation  be- 
tween evaporation  and  the  heat  exchange  at  a 
water  surface.  The  equation  used  here  is 
derived  from  the  relation  found  by  Bowen: 
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*  Values  for  k  are  given  in  Appendix  table  B-3. 


where  Ts  =  sea  temperature  (°  C); 
Ta  =  air  temperature  (°  C); 
and     W     =  wind  speed  (m./sec.). 

Values  of  sensible  heat  exchange  are  gen- 
erally low  in  summer  and  relatively  high  in 
winter  but  nowhere  approach  the  magnitude  of 
heat  flux  through  evaporation. 

Caution  must  be  exercised  in  interpreting 
the  energy  exchange  values  in  regions  having 
limited  observational  coverage.  Small  errors 
in  observation  and  transmission  can  cause 
large  errors  in  some  of  the  computations.  In 
quadrangles  having  few  observations,  consid- 
erable bias  can  be  introduced  by  the  relative 
positions  of  the  reporting  ships  and  their  timing 
with  respect  to  the  calendar  month.  All  com- 
putations presented  below  assume  the  data 
centroid  to  be  at  the  center  of  each  respective 
quadrangle  and  for  the  middle  of  the  calendar 
period  involved.  Energy  exchange  calculations 
are  not  performed  for  5-degree  quadrangles 
having  fewer  than  five  observations  per  month, 
though  summarized  meteorological  data  are 
listed. 

Representative  charts  of  the  seasonal  varia- 
tion of  total  energy  exchange  and  its  components 
are  presented  for  August  1963  and  February 
1964  in  figures  26-37. 

Until  further  experiments  on  the  nature  of 
the  heat  flux  at  the  air-sea  interface  are  com- 
pleted,   equation   models    are    refined,   and  the 


31 


Figure  26. — Incoming  radiation  (cal./cm?/day)  corrected  for  cloud  cover,  August  1963.  All  values  positive. 
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Figure  27. — Incoming  radiation  (cal./cm?/day)  corrected  for  cloud  cover,  February  1964.  All  values  positive. 
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Figure  28.  — Back  radiation  (cal./cm.2/day),  August  1963.  All  values  negative. 
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Figure  29. — Back  radiation  (cal./cmr/day),  February  1964.  AH  values  negative. 
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Figure  30.~Reflected  radiation  (cal./cm. /day),  August  1963.  All  values  negative. 
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Figure  31.  —  Reflected  radiation  (cal./cm.  /day),  February  1964.  All  values  negative. 


37 


Figure  32. — Evaporation  (cal./cm.  /day),  August  1963.  All  values  negative. 
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Figure  33. — Evaporation  (cal./cm.  /day),  February  1964.  All  values  negative. 
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Figure  34. — Conduction  of  sensible  heat  (cal./cm.  /day),  August  1963. 
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Figure  35. — Conduction  of  sensible  heat  (cal./cm?/day),  February  1964. 
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Figure  36. — Total  energy  exchange  (cal./cm.  /day),  August  1963. 
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Figure  37.— Total  energy  exchange  (cal./cm.2/day),  February  1964. 
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accuracy  of  the  synoptic  marine  weather  report 
is  improved,  only  broad  scale  features  of  the 
energy  exchange  will  be  considered.  The  real 
value  of  the  heat  budget  charts  as  presented 
may  be  in  showing  year  to  year  and  month  to 
month  variations  of  each  parameter.  ' 

Until  such  time  as  these  computations  can  be 
improved  more  adequately  to  represent  abso- 
lute values,  they  should  be  considered  only  as 
relative  indices  of  the  magnitude  of  the  total 
energy  flux  which  is  taking  place  at  the  air-sea 
interface.  On  this  basis,  the  data  that  are  being 
produced  routinely  on  a  monthly  schedule 
should  stimulate  further  research  efforts  in 
this  general  subject  area. 
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APPENDIX  A 

Listing  of  the  FORTRAN-62  data  processing  program 
for  synoptic  marine  weather  reports  conforming  to  the 
World  Meteorological  Organization  code  format  FM  21. A., 
International  Surface  Report  from  Ship  in  Full  Form. 
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-COOP > BAA 1-P01 »CL I NE» I /l 0/20/30/0/20/30 .40. 18000 » .NANCY. 

-RINAPY.56. 

(RFLPCOM. 

-FORTRAN.0.1  .56. 

P?JlGf?AM_j^yCfW 

DIMENSION  MFANSFA  (  72.24  )  .  MEANER!  72 »24  )  ."  ISEATEMP  (  10 1 120  )  • 
1 LB_AROM_Q__,  1_2_ 9± ♦__]_«  I_ND_I_R X_t  1_0»_1 2 _0 )  J__J_WI_Npj  ?_Y  '.  1 0  •  1 2  0  )  .  I  W  I  ND VE L  (  1 0  . 

?   12p).  ISFAU0.120) .  IBARI 10.120)  ,  IWIND  (10.120). 
->  TWTNDVd  P. 170)  ,  IRFJT(  10,120)  .  IREJB(1C120  )»  ITE'MPREJ  (  10.120) 

4   .  IBARORFJI 10.120) .  JCLOUD ( 2 . 24 ) »  JDEWPT ( 2  .24  )  .  ^CLO ( 2 » 24 j  . 
5 JDEW(2.24).  JSFA(5.60).  KSEA(2»24)>  KWIND(2.24)  »__K_W_INDVJJ2  .24)  • 

6   KSEATFMP(?.24)  .  ~KW  INDVEU  2  »24  )  .  KW  I  ND  I  RX  <  2~»24  )  »  KWINDI RY< 2 t24  )  . 
7   JSFATEMP(5.60)  »  COUNT  (2.  24).  DIFF  (  2  »24_)_»  KBAROMJ  2  .24  )  .K.BAR  (  2  .24  )  . 

8   XK2(24)»   ALPHA  (2~1  TV   BLANK ( 10 U  Esfl0.40), 

o   ICK10).  IRFJQO).  XK1(?4).  HEAT112.24).  REFL ECT ( 12 >24 ) 

DIMENSION  XREJ(30)  ,XIC2( 10)  »  I  BUFFER  1 ( 1 5 ) . 
J X I C4  ( 2 )  .  ,  I C7 ( 2 ) ,  I C8 ( 2 ) . _  X IC9  <_2  )_»_  1C1.0J.7J  »  _CL0U_D(2Jj>_ 

2   QR(7).  01(2).  OE(?l.  OB ( 2 ) ♦  QH(?>.  QT ( 7 ) .   HEADI120).  AHEAD (360) 
3__j I  C2_(  J  0  ).  BAROM(IO).  I  C_3  (  1  0.)  ..  I  C_4  (  1 0  )  ♦  I.C  5.(  1 0  )  »  I C6  (  1 0  )  » 

4   HEA~D2(60),  AHEAD2I  180)  •  ~HEAD5<24)  •  AHEAD  5  (  72  )  ."XNAME  O")  »  TABLE  (  6  ) 

5     .  AIPTF^P(2) 

COMMON  ISEATEMP,  IBAROM.  IWINDIRX,  IWINDIRY. 
1 IWINDVEL.  ISEA.  I_R_AR_._J_W  I_ND  ,  IWINDV.  I  RE  JT  »__I_REJB.  ITEwPREJ. 

2  IBAROREJ,  JCL0UD~.  JDEWPT.  JCLO. "  JDEW.  "JSEAt  KSEA.'KWIND. 

— JL-JiW  I  NP_Y_»..5_SE.AJ_E_M_P_.__KW_I  N_Dy.EL_.__  KWINDIRX.  KW INDIRY t  _JS_EATEMP_._  COUNT  • 
4   DIFF»lfBAROM»KBAR~ 

COMMON   XREJ.  XIC7>  IBUFFERl.  XIC4.  1C7,  IC8,  X1C9,  IC10,  CLOUD, 

1   QR,  QI,  QE,  QB,  OH,  QT,  HEAD,  AHEAD,  IC2,  BAROM,  IC3,  IC4,  IC5, 
2. TC6»  HFAD2.  AHFAH2.  HFAH5.  AHF Ap5» • XNAMF »  TABLE.  AIRTEMP. 

3  XK7,  ALPHA,  BLANK.  ES»  IC1.  IREJ.  XK1,  HEAT.  REFLECT. 
A MFANSFA.  MFANBAR 

C       **»»*■»»*#»»***»*■******»***»*.*#**#** 


_£ TARI  FS  FORMATS..  . 


1"  F0RMAT(12F6.2) 

_1J_-£0RMAI_L20_A4_1___. 

12  F0RMAT(l?F6.n) 

-+3~FORM*-T-H-AF-8-i-2->- 
\u,    FORMAT(12I6) 

\s,     pnPMAT(ailO) 


C       »♦»♦*#»»»##»*»*»■»**»»##*♦*»»*****»**»*■» 
X JdEAD£.B__CABrL_FnBLMAJL«^.» .....     

16  FORMATdl.    Al.     8X»     II.     Al,     8X,     12,     Al.    7X»     12,     Al,     7X,     12.    8X» 

1 ASj.-1Xj._IAj„6Xi_-2JL1J .. 

r      **««««******»**«******•***•*•********•* 

_C MPHAMFPIf    FORMATS... 

17  F0RMATI6X,     Al,     A?,     2X,     A2 ,    4X,    Al,     2A2,     IX,    A2,    4X ♦    A3.    A2»     13X, 
J J3.LA.lj_2A2.t-_LX±i__3A5J 

18  F0RMAT(4X,  Al,   P(^X,  Al)) 

_C  *  *_iL  ***__******._***_*  *_*__*__*_ *_j!t_JL_^-Jt_?__?-_?__?__*_  *  *****  *.*__* 
C      OUTPUT  FORMATS,  HEADINGS... 

19  FORMATMH1.22X.75HU.S.  BUREAU  OF  COMMERCIAL  FISHERIES     SYNOPTIC 

1WFATHFR  DATA — PACIFIC  OCEAN  /) 

20  FORMAT(1H0.29X.15HLATITUDF  RANGE  . A4»lH-t A4.12X»16HL0NGI TUDE  RANGE 

1  .A4.1H-.A4/  3(1H0   / )  ) 

2a__F0J?Mjii_LiJin^a^9HL0.NiIJ.UD.E_ji?JLj!JLGJ^4.t6.Xi../.9H..LATJ.TUDE_} 

27    F0RMAT(1H0,8X.9HL0NGITUDF,12X.4( A4. 16X) .A4/9H    LATITUDE        ) 
7?     FORMATMH0.SX.9HI  ONGTTUDF.47X.A4.46X.A4     /9H    LATITUDE \ 

C       TFMPS  ONLY... 
?«>  FORMAT(?HQ  A4^  10X»  lO(?X.F?.0»lXt  1H(  t  I2»  1H)  »  IX)    ) 

C      TFMP  RFJFCTS... 
2j6-J0Rl!UlJL2HQ-^4jL-±QAi-J.0J.lAjJi2^Q.iLliLtfJ2.jiliJLXj_F^p.tlXl I 

C      BAROMETER  REJFCTS... 
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2.1-EQ.Bm JAZhD.. AAi..ir_XA.JL0.(.F^O  , 1_X  .E4*0_a  X.LJ 

C  parom,    WIN"    VFCTORS,    AND    '■' I  NO    VEL... 

78    FO    ivaT(?HO    A4.     3X,     3HPAR»     ?X,     10(F6.1tlH(t     12, 1H))»    /       IPX,     2HX-, 

1  3X,     1^(14,     IX,     1H(»     12,     1H),     IX),     /       10X,    2HY-,     3X,     10(14,     IX, 
2 1H(.    12.    1H).    IX),    /      9X»    3HVEL.    3X.    10(14 ,1 X ,    1 H  (  ,    J2,     1H), 

~~ 3          IX)"      )  " 
C  TWO    DfjGRFf    TEMPS. •• 

29  FORMAT (2H0    A4»     13X,       5(     5x~»"F4  .  1 ,  IX  ,"  ~1H l(  ,  ~  12  ,    1H  j",    F5.1)) 

f  FTVF    DFGRFF    FORMAT... . . __ 

30  FORMAT(2H0  A4,28X,8HSEA  TEMP12X,  F4.1,  IX,  1H(,  13,  1H),  40X, 
L F4.I.  IX.  1H(.  13.  1H)/34X»8HAIR  TEMP  12X.F4. 1.1X.1H( , 13. 1H)  , 

2  40X.F4.1  ,1X,1H(  , 13, 1H)  ,/  36X,6HDEW  PT , 1 2X »F4. 1  ,  IX  ,  1H (  » I  3 »1H ) »40X • 
_3_ EAiJLiJ-Xj_lH.(_,_L3jLl_Hl/_ 37X.5HBAROM»10X,F6.1,1X.1H(  . I3.rH)»38X. 

4    F6.1,1X,1H( ,I3,1H)/  40X,  2HX-,  12X,  14,  IX,  1 
5 H(,  T3.  1HI.  40X.  14.  IX,  1H(,  13,  1H),  /  40X,  2HY~.  12X,  14.  1 

6    X,  1H(,  13,  1H),  40X,  14,  IX,  1H(,  13,  1H),  /  39X,  3HVEL,  10X, 
.1 5-6^QjL.lXj-_lJlijL-±3jL-iaj_t_^8A.-_EiL.5-»_J.X.i.AHJ.j..I3-»_A!±Lj__/__3_6X._  6HC 

8L0UDS,12X,F4.2»1X,1H( ,I3»1H) ,40X , F4. 2 » IX , 1H (  ,I3»1H)  ) 

.C TJ.IJ_E.JBAQE.FXlSMAi.jL* . 

3?  F0RMAT(1H1//42X,36H  U.S.  BUREAU  OF  COMMERCIAL  FISHERIES    /, 

] 49X.77H  RTOLOHTfAL  LABORATORY Ll 

?    51X,18H  SAN  D1FGO,  CALIF.      ////, 

_3_ 47X.36H  MARINE  WEATHER  OBSERVATION  SUMMARY _/, 

4    56X,  8H  FOR  THE      /, 

5 53X.14H  PACIFIC  OCEAN //_»_  __ 

6    51  X,  18H  FOR  THE  MONTH  OF    /  ,  53X~,  A8  ,  2H,  ,  I  4",  10  (7  )  , 

7    42X.16H  DATF  OF  RUN...   , A8 , 1 X , I  2  ,2H ,  ,14,   7(/)    ) 

33  F0RMAT(A8,1X,I2,?X, I4,3X,^A8) 

.C D5XGDf._F_QRMAT_.jL. 

35  FORMATII8) 

_3_7__FORMAT(I?A4J         

40  FC^R~MATFiHP,23~X,5H~T ABLE»F6.3»62H — SEA  SURFACE  TEMPERATURE'S  (DEGREES 

1  FAHRENHEIT)  WITH  NUMBER  OF  /  24X ,  56HOBSER  V.AT  I ONS  BY  ONE-DEGREE  OU 

2ADRANGLFS  FOR  THE  MONTH  OF  ,A8,2H.  ,I4,1H.  /) 

_M__F0J?MAT(1H0,26X,5HTABLE,F6.  3.  56Hr_-SE^  SURFACE,  TEMPERATURES  (DEGREES 

1  FAHRENHEIT)  REJECTED  /27X.52HREADINGS  BY  ONE-DEGREE  QUADRANGLE~s"F 

?OR  THE  MONTH  OF ,A8,2H,  ,14, 1H.  /) 

42  FORMAT(lH0»9X.5HTABLE»F6.3t89H~BAROMETR"fc  PRESSURES  (MILLIBARS)." 
IX-  AND  Y-WIND  VFCTORS  (KNOTS),  AND  AVERAGE  WIND  SPEED  /12X,80H(KN0 


2TS)  WITH  NUMBER  OF  OBSERVATIONS  BY  ONE-DEGREE  QUADRANGLES  FOR  THE 

-3-____LLb__QF-----Ai---^JL_-UA-l_i.--_"A 

43  FORMATI1H0.28X.5HTABLE.F6.3.52H — BAROMETRIC  PRESSURES  (MILLIBARS) 
L_UL_IILC_!_-___R-E________L^ 

20F   .A8.2H,  .14, 1H.  /) 

44  FORMATMH0,17X,5HTABLE,F6.3,73H — SEA  SURFACE  TEMPERATURES  (DEGREES 
1  FAHRENHEIT  AND  CENTIGRADE)  AND  NUMBER  /18X.60HOF  OBSERVATIONS  BY 

^IiJQ^DE_GRJEE__OLLA_D_RANC1LEj5_F_OJi_J_HE_  .MONTH.  OF .A8.2H.  .I4.1H.  /) 

45  FORM"AT(1HO,9X,5HTABLE,F6.3.91H--"SEA  SURFACE  TEMPERATU~RES  (DEGRE'eS 
_____NJ__GJ^DEJ._._JLL___l^ 

2INT  TEMPERATURE'S  ("DEGREES  CENT  FgRADE  )  .  BAROMETRIC  P~R~ESSURES  (MILLI 
3PARS).  X-  AND  Y-WIND  VFCTORS  (KNOTS),    /6X  ,  107HAVERAGE  WIND SPE_ 


4EDS    (KNOTS),  CLOUD  COVER  (TENTHS  OF  SKY  COVERED),  AND  HEAT  BUDGE 
J5_I__(£AkQR_I.E_S_/C_E  N  TJ_M  E__E  R  -    /JL0_X  ,  8  &H  SOU  AR E  D  /  DAY).  WITH  NUMBER  OF  OBS 

6ERVATIONS  BY  F I VE-DEGREE  o'uaSraNGLES  FO~R  THE  MONTH  OF    ,A8,2H.  • 
____]________ _/_ . 

46  FORMATI1H1.////.51X.17HTARLE  OF  CONTENTS/9 ( / ) » 1X.74HSEA  SURFACE  TE 
1MPERATURES  (DEGREES  FAHRENHEIT)  WITH  NUMBER  OF  OBSERVATIONS   //1X. 


243HBY  ONE-DEGREE  QUADRANGLES  FOR  THE  MONTH  OF   .A8.2H.  .I4.1H..40X 
--2__U-iIAe]^ES--__-_P_l-_f.6____zVy__^ 

4»  X-  AND  Y-WIND  VECTORS  (KNOTS').  AND  AVERAGE  WIND  SPEED  (KNOTS).// 

51X.71HWITH  NUMBER  OF„OB_SE£.VADONS  J_Y  ONE-QEGREE  _QUApRA_NG_LES.F_OR__TH_ 


6F  MONTH  OF   ,A8,2H»  .I4.1H.  12X  ,  13HT ABLES  2 . 001 -.F6.3 . // //// IX, 72H 
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7SEA  SURFACF  TEMPERATURES  (DEGREES  FAHRENHE I TJ  AND  _NUMBER__OF__OBSERV 
8ATI0NS   77TJTi43HBY  TWO-DEGRFF  QUADRANGLES"  FOR  THE  MONTH  OF   ~»A8> 

92H,  .14, 1H.   40X,1 3HTABLES  3.001-   F6.3/////) 

47  F0RMAT(45X,24HNUMRER  OF  OBSERVATIONS   I  5 /45X .24HNUMBER  OF  REJECTIO 

INS I5M5X.24HSUMMARIZED  OBSERVATIONS  15//////// 

2    42X.33H* PROGRAMMED  BY. 7-  MARVIN "w."  CLINE   // 

. 3_ 4J_X,__l_8H._PROC.ESSE_D_  BY.... 3A8) 

4P  FORMAT)  14X.  2( 

1  MX,  39HQ(T)  =  Q(I)  +  Q(R)  +  Q(B)  +  Q(E)  +  Q(H)  )<  /  14X.2  (  11X. 

2  F4.0.  3X.  F4.0,  3X,  F4.0,  3X,  F4.0,  3X.  F4.0.  3X»  F4.C)    ) 

50  FORMAT ( IX »85HSEA  SURFACE  TEMPERATURES  (DEGREES  CENTIGRADE)! .AIR._TE. 
1MPERATURFS  (DEGRFFS  CENTIGRADE).  //1X.84HDEW  POINT  TEMPERATURES  (D 

2FKRFES  CENTIGRADE.  BAROMETRIC  PRESSURES  (MILLIBARS)*  X-  AND  //1X» 

383HY-WIND  VECTORS  (KNOTS  )  ,  "a've'RAGE  Vl  ND    SPEEDS    (KNOTS).  CLOUD 
4C0VER  (TENTHS  OF  SKY  / / IX , 76HCOVERED ) .  AND  HEAT  BUDGET  (CALORIES/C 


5ENTIMETER-SQUARED/DAY) .  WITH  NUMBER  OF   // IX .  57HOBSERVAT I ONS  BY  FI 

fLYEjrDi:CilEL-ilUADR-A.NJiL-E.S-.FP_R__IH_E_jMOMJH.OF __________  . I4.1H.  . 

7    26X.13HTABLFS  4. "01-  F6.3/1H1) 
51  FORMAT(1HO»29X.15HLATITUDE  RANGE  »A4t 1H-.A4»12X» 16HLQNGITUDE  RANGE 
1  ,A4,lH-,A<t/  /) 


CONSTANTS... 
A~15-_C 

n=15.  n 

J2Q_6?...I.=lt6. 

f,o  TAn|_F(  I  )  =  I 
'K  =  10 


JL  =  20 
|M_=3_o 

JJJJ=o 
__-__--_=__ 

I"»  =  10»*3 

14=10**4 


T  S= 1 0**^ 

____________ 

I7=10#»7 

_L___2__?____. 

JQ  =  1  n#»«? 

110=10**10 


111=10**11 

L______Q___I2 

lli=To**i3~ 
c     *♦♦#»**■»*#»***»****************< 

C      RFAD  IN  PERM.  TABLES... 

RFAD  10.  XK1  .  XK? 


RCAD  11.   HFAO,  AHEAD.  HF*n2»  AHFAD2 
RFAD  37»  HFAD5.  AHFAD5 

RFAD  12,  HFAT.  RFFLECT 

________13.__JE__ 

C   ##**»»»»*♦»#»»»♦♦»**«****»*********** 

£ RTAD  HFADFR  CARD  AND  S^T  IIP  LIMITS... 

C       PULL  PRECEDING  CARD    REPLACE  WITHOUT  STATEMENT  NUMBER 

6969  READ  INPUT  TAPE  10«16»  LIMLATN. DIRN_t__L  I MLAJS » D  IRS  »  L  I  MLONGW  • 

1    DIRW.  LIMLONGE.   DIRE.  NO.AMONTH.  IYEAR.  MARKT,  MARKB 

READ  INPUT  TAPE  10.  33»  XDATE.IDAY*  JYEAR.  XNAME 

C      RFAD  MONTHLY  TABLES... 

RFAH  INPUT  TAPE  10.14.  MFANSFA 

RFAD  INPUT  TAPE  10,15.  MFANRAR 

lF--iJ_JllLAm-6-U4.DH0jLin.Qil.UilO.fi. 

1000  IF  (DIRN-1HN)  1002.1001 ,1002 

.-4-0ill~LJ--_AIJ--A.Q-XJ--U-I-J.-4J- -— 

GO  TO  1007' 
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_Lan2_IJL-JLOiRtirJJdS.L_J.Qil4-t.LaQijJ_on.4i. 

1003  L!MLATN=60+LIMLATN*in 
GO  TO  1QP7 


1004  IF  (DIRN-1H  )  1006,1005 »1006 
.iQi25._jj..xiJL_L._Jj-J.-__QD---_LQ_ll_'_lP_Q6_ 

1006  LIMLATN-0 
_iQ_27-_I5__±LJL(____J_Sr__J__l±07.iJL10.7jj.0J.3. 

Lf°7    IF    (0IRS-1HS)     loog, 1008,10P9 

1PPR     I   TMI  ATS-6P+I   TM|  ATS»1P 

GO    TO    1^14 

__Li'_?__J-F__LQ_L_S.rJ___Ll.-IC_Ll__l_110_j_._0_._l 

lOl"  LTMLATS  =  6r,-LIMLATS*10 
GO  TO  1014 

101]     IF    (0IRS-1H    )     1013.10~12.1013 

1Q12    IF    (LIMLATS)     1013,1010,1013 


1013  LIMLATS=12P 

1014  IF  (LIMLATS-LIMLATN)  1015,1015,1016 

1015  LIMLATN-0 

LIMLATS=120 

1016  IF    (LlMLONGW)     ln2?  1 1 1 16  .1115 
111^    IF     (LIMLONGW-18)     1116,1116,1022 


1116     IF     (DIRW-1HW)     loi p, 1017,1018 
1017    l.TMLONGW=LI^LON(_W»10 

GO    TO    102  3 
101P    IF    (DIRW-1HE)     1020,1019,1020 
1019    LI MLONGW-360-LI~MLONGW»10 

GO    TO    107  3 


10?0     IF    (DIRW-1H    )     1022.1021,1022 
.iQ_____J_L__LUL-!__QN_i_-rl_il__l__22.ili___lj_LQ_L2_ 

102?    LIML0MGW-70 
.i02.3__J5_-_a±-_jQN_G___l_8_l__l__2J.j.LI2i»_LQ.2.9. 

11?^      IF     (DIRF-1HF)      102^,1024,1025 
1074     I   TMI  QNGF  =  360-I    1MIOMGF-10 


GO  TO  1O30 

.i___5__i.F_-__D.iJ__ir___yi_i_z__aa2-__--i2-i 

1026  LIMLONGF=LIMLONGE*10 

_____Q_„r„0_ 

1027  IF  (0IRF-1M  )  1029,1028,1029 

1028  IF  (LIMLONGF-18)  1029,1026,1029 


1029  LIMLONGE=180 

1030  IF  (LIMLONGW-LIMLONGE)  1032,1031,1031 

1031  LIMLONGW-70 

LJ„LjO_G_=l_0 

103  7    LONGCHK=LIMLONGW+10 
I   ATM1-I    TMI   flTM-H 


LATS1=LIMLATS 
-U_I____LA__L_._L___ 

LATS2=(LATSl+l)/2 
-LAI-i-_=_iL_I____t__L__i- 

LATS5=(LATSl+4)/5 

IF  (NO)  1033,1033,2033 


2033  IF  (NO-12)  1034,1034,1033 

1033  N0=1? 

1034  CONTINUE 

_ro____L=__ 

ASSIGN  6000  TO  NOT 
AOOO  TF  MINTT.   1.1  AQQP.ftOPI  .7003.7000 


60O1  BUFFER  IN  t JK , 0 )  ( I  PUFFER!  ( 1 )  ,  I BUFFER1 ( 1 5 )  ) 

C_ •_»_»._»_»....»  »_»  »  »,*  *  *  »  *  *  *  *  *  *  *  »  .«  *_.*..*.  *_*  *..*..*  .*.*—*--  * 

r      ZFRO  ARRAYS... 

-_l__3i__AL____-Efii_-_L_-_.9__Q_L__-SEAT_L_'P_l 

ASSIGN  6002  TO  NOT 
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C  15fY6^~1)mTTnT6~7lP^A~nT^ 

1^3     irouNTs" 

60^2     IF     (UNIT.     JK )     6002,     104,7003,7000 
104    OFCOPE    (    80,    17 ,     IPUrTERlJ ALPHA.. 

of^odf   (   6a  »  iT» ~T bTJf'fer'i  j   r"l ank  " 

lO.^.BU F_FF_R__T  N___(_JK_,  0J___(_ Jjiy  FJfFRjJ  1  )  ♦  I  PUFFERl ( 1 5 )) 

C      TESTS  FOR  XS  AND  BLANKS   ALPHA  AND  NUMERIC 

106  IF  (ALPHAM9)-5H»»«*»)  107.204,107 

107  1=0 

_il£L_L*l±l . 

IF  (BLANK(I)-IH  )  103.108,103 

108  IF  (1-10)  5107»5108.5108 

"5108  IF  (D"OLT(  ALPHA!  1  )  ,1)  V   103  .  1 10"t  103 

Un  DFCODF  (8,  35,  ALPHA  111)   IOCTANT 


I0CTANT=I0CTANT/I7 

lE.AHQLIJALPi±MZl.\.Z.U..l9.1t.U2.i±^. 

112    "cCODF     (8,    35.     ALPHM2))        ILAT 

LLAIr-LULLZlfi 

IF  (DOLTt ALPHAI3 ) .2) )  103.114.103 
114  DFCODF  (8,  35,  ALPHA(3)1   I  LONG 


IL0NG=IL0NG/I6 
C   *»»*#»»»***»*#♦**»#*♦»»#*»***♦**»#***** 

C  TOMPUTXfTOfTTJFniONGHATIirTTf^  "  5~  d¥g" 

C  ARRAYS... 

117  IDUMMY=IOCTANT+i 

GO    TO     (118,118,118,118,103.119.119.119,119.103).     IDUMMY 

118  L4T=60-ILAT 

.QfULQ-AJDiL 

110    LAT=61+ILAT 

$10    IF    (LAT-LIMLATN)     103,103,601 

6r\     IF     (LAT-LIMLATS)     122.103.103 

12?     M|   AT=(I   AT+4W5 

*LAT=(LAT+1 ) /2 
.„_L22_JaCLjrJLJLl25.JLl2&.*J^.2*.m*^^ 

126  IF  (ILONG-80)  127.103.125 
. 12-T.J.UlNJirJLLO.NG.tl.oC. 

125    LONG=10-ILONG+ILONG/10»10 
ti  nMG=Ti  nwr,*) 


124    ML0MG=(L0NG+4) /5 
ja.Qfc!£=-LL.0Airvfel-)-/-2. 

GO  TO  135 
_U3_I£_l-LLflNG=anj_ lJ.nj.lH3jJ.3_L- 
130  ILONG=ILONG+100 
T»1  |  nMr,=  n-T|  nMr,-T|  nwr,  /m*i  o 


ILONG=360-ILONG 

.go— 10.-12.4. 

135  jlong=(  il0ng..4  )/5 

C !L_JL_«_JLJLJI_JtJt_JLJt_Jt_A_JLJL_*_S_JLJLJS_*_ft_JLJL_±_i 


C      IS  CARD  IN  LIMITS... 

TF     fMFANSFAf.il  ONG. Ml  ATM     103.999.999 


090  IF  (LONGCHK-ILONG)  103.136.136 

. ia&^F_iiLaMrL-=iJ3N£c±tKjLir-i--inajjLP3_i-i32. 

c    »»*»»♦♦»»»»♦»»#**♦»♦*»#*♦**»****»**»#** 

C ZFfiO-I.MDlCATaR— SEARCH.^..*- 

137    T=7 
5137    1  =  1  +  3 


IF  ( ALPHAf I 1-1H0)  138.139.138 

. — l3^-;F-I4-=ri64--5-1.3Jlji73.jl75- 

C   »#»♦*■*■»■»»***«■***»»****»*»*»»*******»*** 

C SEA— TEMe£fiA.TURE_.SHBRQUIiME*..j. 

130  IF  (DOLTt ALPHAI9) ,2) )  175,141.175 
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lAl__D£.01DfL-lfljL_3L5i_-Ai.eH&JL2J-l HEMP 

ITF^P=ITEmp/I6 
IF     (DOLT( ALPHA! 1  +  1 ) .2)  )     175.142,175 


14?    DFCODF     (8»    35.    ALPHAII+11)        IDIFF 

U1JJF^=IDI  FT/16 

■  SIIP-ZFRO"  ANO    OTHER    JUNK    SUBROUTINE. 

___I_F___(_LO_NGCHK-1801  9550  .  9550  .9500 

9  50  0  IF  tNO"-3~)  "95~10»9~5"l07l43""~ 

951"  IF  (LAT-20)  9530 .9530 .9520 


9520  IF  (LAT-101)  143.9530.9530 
Jl53_0..JJL_XIJ.EaPr.5D_L.J.^3.»J.^A»9.54_0_ 

9540  ITEMP=50-ITEMP 
G.O.T_0__14_3 

955  0  IF  (NO-3)~  9 560, 9 560", 80oi 
956"  TF  MAT-?")  958" . 9580 , 9570 


957"     IF     (LAT-101)     8001  .9580, P5fl0 
J25J10__jJL_jLIJLE&PrJ5£Ll-JioJHjL8ILGki9A9J2_ 

9591-1     IF     (ITFMP-79)     9600.9600.9610 
.-9_6_C0.-J-tEM-e=_50--rHEMP_ 

GO    TO    143 

9610     TTFMP=TTFMP-1"0 


GO    TO    147 

_aojll„IiL.XCJIlr3-l-^7-L0j.9.11£.i5.ZDJ2 

9700  IF  (NO-12)  9730,9710.9730 

-J5Ii£,--J£._±LAIrj=tr.L__97-2.QjL?J-2.Q.»  i.73.0. 
972"  IF  UTFMP-3")  143_, 147,147 
973"  TF  (ITFMP-36)  143.147,147 


143  IF  (IDIFF-50)  145,144,144 
14_4__jj_EM.PAlJ_EMPi!J.DiiJLDJ_Ffjr50J.*l_Q./2. 

XDUMMY=IDIFF 

PJF_FJM.Lpj^G_»MJ_AJ_l=J}JFFJ_MJ_^^^ 

GO  TO  146 

145    JTEMP=ITEMP»lP-IDIFF»10/2 


XOIIMMY=IOIFF 
...PJ_LFi!iLPJ^G_tMJ.AJ_i=_DJ_FF_(_MJ.g_NGj_ML_AT^r_X 

STATE^fTTt'T^Y'c'oNVERTS    TEMP    FROM    C    TO    F . . . 

AAt..jyJE^PiLUlE}l'P*A3J2X/}.rl£lA3.Z0 

GO"   TO    15" 
147     IF     (IPIFF-50)     148.149,149 


148  JTEMP=( ITEMP-IDIFF)»10 

X.nUMMY.=  I_QiEF_ 

0IFF(MLONG,MLAT)=DIFF(MLONG»MLAT)-XDUMMY/1.8 
GO_TO..150 

149  JTEMP=(  ITEMP+IDIF"F-50)"*10 
XDIIMMY=IDTFF 


DTFF(MLONG.MLAT)=DIFF(MLONG.MLAT )+( XDUMMY-50,0) /1.8 

15_Q__cp_UMJ_MkQM.ilLAlX=^QyJlIJ-Mi.QiiGjLMkALL±_l*0. 

c  compare  with  mean  temp... 

!l=LAI-LAI/_LD.*l_n_ 

IF  (N)  103,151.152 
T"  N=10 


152  IFINK=MEANSEA< JLONG.MLAT) 

LiMJLI.=  l£.lNJl/JLCilO 

TF  (IFINK)  103,160.153 
_113_Jf_.JLN_-3J__JJi4_U6_ljJ.62 

154  NN=3 

155  NNN=-1 


157  IF    (MLAT-1)       103.161.158 

-i&3-jm=a 

NNN^l 
-156— IE.  (Mt.AT_-2.4J__1.5_a.iLlfajJ.P.3_ 

158  IF  (MEANSEAt JLONG.MLAT+NNN ) )  161.161.5158 
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5158    IMEANgIFINK-LIMIT»1000-(IFINK-LIMIT*1000-MEANSEA( JLONG.MLAT+NNN )+ 
1      MEANSEAULONG^MLAT+NNN")"7rob6»10od )*XABSF  (  NN-N  )  /5 

159  IF  (XARSF(JTFMP-IMFAN)-LIMIT«10)  160,160,167 

160  I SF A  TEMPI LONG, LAT ) = JTEMP+ I 5E A  TEMP ( LONG tl AT ) 

ISEA(LONG.LAT)=ISEA(LONG»LAT  L+.1 

JSEATEMP(KLONG,KLAT)=JSEATEMP(KLONG,KLAT)+JTEMP 

JlSE&.«.KL.QN.GjJiLATJ.rJIS£.6I.K^QRQj.KL.ALl+A 

KSFATEMP(MLONG,MLAT)=KSEATEMP(MLONG»MLAT )+< JTEMP-320 )* 100/180 

t^FflfMI  ONfi.VI   AT1±ICSFA(MI  QMG.MI  AT  1+1 

HO  TO  170 

...±6±.JL^I.&n=l£l!ilL-Lmn*.lSlQ0. 

GO  m  159 

..._L£i2_.LF..iJlrAl-_16.4jJ.6-LjJL6_l 

164  IF  (N-6)  165, 166,166 
1f>5  NN  =  ^> 


NMN=1 

GO  TO  156 

1*66  NN  =  8 
NNN=-1 

GO  TO  158 

RFJFCT  ARRAY  PILL... 


167  IF  (JTFMP)  1 ii ,9167,9167 

_9_L&.7__LF--1-»1&RKJ_1__UA7j.1-13j.1.L6.J- 

1167  DO  168  J=1.3 
IF  tIREJT(LONG«LATl!-J)  169«168»16B 

168  CONTINUE 

I.JJ.J  =  J.JJJt] 


GO  TO  1"3 

169  ITEMPREJ( LONG* LAT l=ITEMPRFJ(LONG»LAT)+JTEMP/10»10««t (J-l)«31 

IREJTI  LONG, LAT  )  =  IREJT  (  LONG»'LAT  1  +  1 
-■1J--IJJLJJ-.JJ+1. 

GO  TO  103 
r  ************************************ 


C      DFW  POINT  SUBROUTINE... 

170  IF  t  DOLT  (ALPHA  (1+2 1.2)  ) 175»1?_1.»  175. 

171  DECODE  (8,  35,  ALPHAII+2))   IDEWPT 
IDEWPT=IDEWPT/I5 

8002  IF  ( IDFWPT-990)5170»103,in3 

517^  IF  (ITEMP-36)  5171,5172,5172 


5171  IDEWPT=IDEWPT*180/100+320 

5172  IF  (XABSFI IDEWPT-JTEMP 1-200 )  172*172*175 

172  JDEWPT  (  MLONG»MLAT  )  = JDEWPT  (  MLOiMG  »MLAT  )  +  IDEWPT 

JDEW(ML0NG*MLAT)=JDEW(ML0NG*MLAT)+1 

j-  ************************************* 

C CLOUD  COVER  SUFRQUTINEtf 


1"T5  IF  (DOLT(  ALPHA(A)  ,1  )  )  181,177,181 

lI.7-X!EnOIlE-.lfl.»._35j.-ALP.bAL4LLJ- LCL.QUD— 

ICL0UD=ICL0UD/I7 
..BJ3£L3.__lF_iJ.a_011Dr.9.J„La.Oj.lJ-ajLlfi_L 

178  IF  (DOLT( ALPHA(7) ,2) )  181,800,181 
800  DFCOnF  (8.  ?5.  AI.PHA17)  1 I  V  I  S, 


IVIS=IVIS/I6 
IF  (IVIS-90)  2000»2001»179 

179  IF  (IVIS-93)  2001,2001,181 

.2-on.Q_.iF—  ij.yis-.5j ^on.ij.2iiQj_iiJai 

2001     ICL0UD=8 
1  flO    .  lc-1  Oliri  t  Ml  ONG.MI  ATl=.JrL0ur)(ML0NGtMLATl+l  T  CLOUD*!  00/8  ) 


JCLO(MLONG,MLAT) = JCLOt MLONG»MLAT )+l 

£ JL.JL_^-^_*_-*-.*-.*_Jt..*_Jt-A-it.JL-*-_#--*--*--ft.-*^^ 

C  RAROMFTFR    SUBROUTINE... 

— t5.}.__|.f__^JiOL-I-l-A-L42HA4-aj-^34-i— 103.*J,83.«J,Oi 

183  DECODE  (8,  35,  ALPHAI81)   IPRESS 
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_J5.R£.S.S=J_P-R£SS/-_5 

8004  IF  ( IPRESS-500)  184.103.1P5 
184  TPRF.SfisTPRFSS  +  IPPOO 


GO  TO  186 
__La5__lP_B£Si=lER£S5±£D-0-0__. 

186     IFir'K=MFANPAR(  JLONG.MLAT) 
LL«LLI_=If_LNJC/10_1D_o 

TF     (IFINK)     10?, 189. 6500 
6500    N=LAT-LAT/10»10 


TF     (N)     103, 7™. 711 

.ID.rL_M3.lG. 

701     IF     (N-3)     702.6700.704 

.7J12_.M=-1 

NNN=-1 

GO    TO    6501 


704    IF     (N-8)     706.67^0.705 
.7.05— M=-8 

MMN]  =  1 

-GO.T_0__7P.3 

706  IF     (N-6)    "707,7^8.708 

707  NM=3 


NMN-1 
GO    TO    187 

7"8    NN  =  8~ 

NNN—1 

GO    TO    187 
7^3     IF     (MLAT-24)     187,6700,103 


6501     IF     (MLAT-1)     ln3.6700,187 

6700    IMFAN=IFINK-LIMIT*10000 

GO    TO     188 
187    TF    <MFANPAR( JLONG»M.AT+NNN) )    6700t6700»5187 

5187     IMFAN~=IFINK-LIMIT*10000-UFINK-LIMIT*10000-MEANBAR<  JLON~G  .MLAT  +  NNN  ) 
1     ^.VIFANqARtJI  ONn.  Ml.  AT+NNN  1/10000*1  0000  ) »X APSF ( NN-N ) /5 . 


188     TF     (XAPSFt IPRFSS-IMEAN)-LIMIT*10)     189,189,190 
.li9..1ilfiilQWLL.Q^inJi-AJLijLlP.4RO^LLJlNjQikATaiLLPJiF.aS 

IBAR(L0NG,LAT)=IPAR(L0NG.LAT)+1 
|^AiLQM_lML0J^G_tM^AT_!5_KB_ARJD^^^ 

KPAR (MLONG,"LAT)=KPAR(MLOMG.MLAT )+l 
GO    TO    193 


C       RPJFCT  ARRAY  FILL... 

1.9.0„I5_-±yASl<L,5J..lJ.90.,1.03_1i_L9_Q 

1 19"  00  191  J=1 ,? 

lZ_JLI_R.EOP.J_L0_MGjLLAIJ.-J-L_J.2J>.iJL2JLil.,?l_ 

191  CONTINIIF 
IJ.I.I  =  „I„I  .1,1+1 


GO  TO  103 
-_:i5.2-lH6ilQREXLlJ2Nnj.ULLl=J.BAR£Rj:JJLL0 

IREJB(L0NG,LAT)=IREJB(L0NG.LAT)+1 

LUJ.r_LJJJ+-l 

GO  TO  103 

C      WIND  SUBROUTINE. ■• 

15.3— JJL_XDIU.J-LAJ_PHAX6X»ZIJ— 1113^.15.5.  jl10J3 

195  DFCODF  (8,  35,  ALPHAI6))   IVEL 

4-*£E-_-=-i-UF4,V-T-6 

196  TF  <OOLT( ALPHAI51 .?) )  103,198.103 

198  nFCODE  lfl»  ~i5t  ALPHA15)) 1  P  I  R    


iniR=IDIR/I6 

JX_jLIIlIR.l-J.£l3.iJ.9i»±.99_ 

197  IVEL=0 
G0-TO-2O.2 

199  IF  (IDIR-51)  202.20D.20n 
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2  °  °    IF    (ID I R  -  9  9  1    2  01, 103.103 
TO  1    !  n  I R  =  I b  I  R~-  5  6~~ 

TVFL=IVFL+lon 


?n7     TriIt?=iniR*10 

V/_L=  IVFL 

x^l  imm  y=  iTTir 

IX  =  SINF(Xr)UMMY»0.0174)»VEL*(-l-.Q) 

I Y  =  COS F  ( XDU MM" Y » 0 .  "o  1 7 4~)  * V EL*1  - 1 .  0  )  "" 
KWINDIRX(MLONG.MLAT)=KWINDIRX(MLONG.MLAT)+IX 
KWINDIRY(MLONG.MLAT) =  KW I Nn I RY ( MLONG .ML AT )  ♦  I  Y 

KWIND(ML0NG»MLAT)=KWIND(ML0NG»MLAT)+1 

IWINdIr~X(LONG»LAT  )=IWINDIRX(LONG»LAT)+IX 

IWINDIRY(LONGtLAT)gIWINDIRY(LONG»LAT)+IY 

IWIMD(LONG,LAT )  =  I W I  NO < LONG, LAT > +  1 
7^  IWINDVFL(LONG,LAT)=IWINDVFL(LONG,LAT)+IVEL 


H.'IM^V(LOMG,LAT)  =  TWINnV(LOMG»LAT)  +  l 

KW I NPVFL  (  MLONG iMLATJLgK-WJ^PVFJ_JJ^lJ3NG iML AT  1+ IVE.L 

KWINDV(ML0NG,MLAT)=KWINDV(ML0NG,MLAT>+1 

_6Q_IQ.JLQ3. 

»*#»»»*♦♦»#■*»*♦*»*****»**»■**»■»#»*»»♦** 

START  OF  HFAT  EQUATIONS  AND  OUTPUT  SUBROUTINES 


C      ONF  DEGRFF  OUTPUT  (TFMPS  ONLY) 

C NOTE...  TWO  REELS  ARE  REQUIRED  FOR  OUTPUT  IF  IT  IS  INDICATED  ON  THE 

C    ~  HFADER  CARD  THAT  REJECTS  ARE  DESIRED. 

204  I  F  _  ALPHA_(_7_0_  -5H#***»  )  5204  »  5 24 7 _»___2_0 4 

"~52  0  4  IF  (~(  T  ABL  Ell)-  lT6  )  *  3  9  .  O+TTABL  E  <2)-2.0)*57.0+T  TAB  L  ~E  (  3  )"-  3  .  0  )  *  3  V.  0  + 

1   (TAPLF(4)-4.0)«71.0-A)  5706,5205.5205 


«705  CALL  CHANGE  (JL,A) 

57106  IF  (  (TABLE(5)-5.0)»39.0+(TABLE(6)-6. 0)*39. 0-B)  520.8t5_207>  5207 
5707  CALL  CHANGF  (JM,R) 
_5J._L8__M__LD_NJ___H._-9_ 

TARLEI  1  )  =  TARLFU  ) +0.001 

TARLF(5)=TAPLF(5)+O.OCl 


L1NF=0 
AQ0Q_j!iRlT£_P,UJ_PUI-_I&P.EL-Jki..l2.t 

WRITE  OUTPUT  TAPE  JL,  40,  TABLE! 1),  AMONTH.  IYEAR 
WR_LIJI.0UT.PUI. _t_J_E__JL_»__2C_____EAD___AIN  1J_»_  H_AD____A_TAlJ_t_  AHEA.D  (M  )  j, 

1   4HFAH(L0NGCHK) 
WRTTF  OUTPUT  TAPF  JL.  71.  (AHEAD(I).  I=M,LONGCHK) 


4nr>7  TF  (MARKT)  70^,706,705 
.-Jl!:ji_w?_LIf-nUJ-PlLI.-t4.P.F__J_M_i__19.1. 

WRITE  OUTPUT  TAPE  JM,  41,  TABLEI5),  AMONTH,  IYEAR 
WRTTF  OUTPUT  TAPF  J^t    2Q»  HFAD(LATN1>»  HEAD(LATS1)»  AHEAD ( M ), 

1    AMFAri(LONGCHK) 
WRTTF  OUTPUT  TAPF  JM.  21,  (AHEAD)  I  ),  I=M,LONGCHK) 


?06  HO  714  J=LATN1,  LATS1 
_QO_.2J-2„IiLljLlO_ 

IC1 (  I  )  =  1SFA(  I ,J) 
_LF._1JLCj_1J_LJ-_2j_0jl2J_Qjl2P_7 

?r,7  T<-?(  I  )  =  isfatfmP(  1  ,j)  /IC1  (  I  ) 

TF  ( TC7( T )-TC?( I 1 /10»10-5)  708,709,709 


7"P  XTC7I  I  )  =  IC7( I ) /in 
GD__T0_5____n 

7^9  XTC7(  I)  =  (  IC7( I  )+5) /in 
J30„Tn__100_Q 

710    XTC7I  I  )=0.0 
SOOO     TF     IMtPlTTl     7TT.717.7TT 


711  IRFJ( I*^-?)=ITFMPREJ< I »J)/I6 

IELEJ.LJ«3---lX=J_LIJIMP_RJlJ__LJLjJ---JJiFaiJL*3-r2_L*_-lA___iy-l3_ 

IRFJI I*3)=ITEMPREJ( I.J)-IREJ( 1*3-2)*  16   -IREJt 1*3-1 )* 13 
XEJ£J-LJ-*.3^--2J-=iEEJ-i-I-*-3.=-2-l 

XRFJf 1*3-1 )  =  IREJ( 1*3-1  ) 


^ 


.XSLE-ll±*3-l=-l9JiJL.Ll*.3±- 


7i  7 


rnMT  \  mi  ic 
ITNC=LINF+1 


9000 


30)  400^,4005.9000 

TARLFtl )+0.001 


IF  (LINE 

LTMF="l 

9  00 1  wd  I TE  OUT  PUT  TAP E __J L  »  19  . 

'  "write  output  tap"f  JL,  46, 

WRITE  OUTPUT  TAPF  JL,  20, 


TARLE(l).  AMONTH,  IYEAR 
HEAD(LATNl),  HEAD(LATSl),  AHEAD(M) • 


1   ahc an (LONGCHK  ) 
WRfTf  OUTPUT  TAPE  JL»  21»  (AHEADtI)»  I=M»LONGCHK) 

DdM  IF  (MARKT)  9004,4005,9004 

Jli?Jl&„TABJ_Fj.5Jj;TjAf^fj_5j_+^tjinj 

9niq  WPTTF  OUTPUT  TAPF  JM ,  19. 

MPTTF  OUTPUT  TAPF  JM  ,  41,  TARLFI5),  AMON  T  H  »  I  Y  E  A  R 

HFADI LATS1 ) , 


wpitf  output  tape  jm,  20,  headilatnd, 

L    _AHFAUlLO_MG_CHj>r  j 

write  output  tape  jm,  21*  (  aheadh  i  »i=m»  longchk i 

4005  WRITE  OUTPUT  TAPE  JL,25» HEA_DLJ1j»___(_X_I.C2  (  I  )  »  IC1  <  I  )  » 

4007  jr  (Mari^T)  213,214,213 
71?  WQTTF  OUTPUT  TAPF  JM,26,   HEAD(J), XRE J 


AHEAD(M) , 


1=1.10) 


714  CONTINUE 

_*  *  *  *  _*_*__5  *._*_  *_  * *  _*_J 

PFMAJMnrp  OF  OMF  ncG9FF 
J  o^L rX? )- TA ? L  F(7)  +  n.O01 

TAQLF(A)=T^OLF(6)+^.Onl 

1  T"c  =  o 


_  *  *  *  #  -I 

OUTPUT... 


jLIMF=n 

22 1*  WP. ITE  OUTPUT  T APE  JL»  19» 

WRITE  OUTPUT  TAPE  JL.  42.  TARLEI2).  AMONTH,  IYEAR 

WRITE  OUTPUT  TAPE  JL.  20>  HEAD1LATN1 ) .  HEAD ( LATS1 ) .  AHEAD (M ) . 
1   AHEADfLONGCHK ) 

WPITF  OUTPUT  TAPE  JL .  21.  (AHEAD(I).  I=M.  LONGCHK) 


471  «i  IF  (M4RVP)  715,716,215 
71*  WPITF  OUTPUT  TAPF  JM ,  19 

"  ~Wf»IT?  OUTPUT  TAPF~JM,  43 

w5.IJLF__Qy_IpJJL_TAPf._  J_^_»-.  20  __ 

1   «MF/<n  f  l  OM^THie  ) 
WRITF  OUTPUT  TAPF  JM.  71 ,  (AHEAD!  IJ»  I=M»  LONGCHK) 


TABLEI6),  AMONTH.  IYEAR 
HEAD(LATNl).  HEAD(LATSl).  AHEAD(M). 


716 


no  77«,  J=L.ATM1  ,  LATM 

00  ??4  1  =  1.10   

!<"?(  I  )  =  TPAR(  I, J) 
IF  (IC7  UIJ  2  1  R  .  2  1  R  .  2  1  7 

71 7  RAROM( I )  =  IRAROM(  I  ,J)  / rC2(  I  ) 

RAROMI  I  )=PARQM(  I  ) /I  0.0 

GO  TO  500] 

__?_\.5._?A?^J_IJ_=2.«._L— 
s^i  t^ (i ) =iw i"NrT(  f  ,~j i 

TF  (  IC*m  1  7  7"  ,  7?_o  .  ?A?__. 

7io  rr«i(n  =  TWIMOTRX(T,j)/lC3(T) 

TQ4(n  =  TWlNOTRY(T,J)/ir?(T) 

r,0  TO  "5  on  7 

??n  Tr4(l)=o 

Tr5( I )=o 
500?  I("6(I  )  =  TWINHV(  I ,  J  ) 

IF  (ir6(Il-r  222.222. 221 
221   TWINDVFLt  I  .J)  =  IWINDVEL 1 

no  to  son? 
—22.2-1  HUMDVELil-t.  J.)  =_0___ 

S""1  Tc  ("ARKR)  273,774,223 
_  -25a__T_R£_H-L*2--l  )  =  I  SA.5QBE  Ji  I  ,  J  )  /  I  6  ... 

IRFJ(  1*7)  =IBAPORFJ  (  I.J)-IRFJ{  1*2- 


.J) /IC6(  I 


11*16 


58 


XPFj( I»2)=TRFJ(I»2) 
774  CONTINUE 


LINE  =LINE+1 
_lLJ_N£=_iU_Nf.+] 

IF     (LINE-12)     9010,9010.9007 

9007    TflRLEC?)  =  TARLF(?)+0«001 

LINP=1 

Pfffl  WPTTF  OUTPUT  TAP.F  Jl  .  19. 


WITF  OUTPUT  TAPE  JL.  42.  TABLE! 2).  AMONTH.  IYEAR 

_W!UI-l_P_yTJ_m__lA__E___LL_____a_^^^ 

1    AHFAniLONCCHK ) 
WRITE  OUTPUT  _TA_PF__J_L.__2_l_.__(_A_HE_Ap(  IJ.  I=M,_  L0NGCHKJ_ 
"o ni7^"fF~7 mTrTr  r~9 O  11,4011.9011 

9011  IF  (JLINE-30)  4011.4011.9012 

901?  JLIMF=1 


TARL F(6)=TABLF(6)+0.001 

WP  f f~E "output  ~f  apT'Tm",  19, 

WRITE  OUTPUTTAPF  _JM^,__43,  TABLE  (_6  )_._  AMONJH.  IYEAR 

WRITE  OOTPUT  TAPE  JM»  20,"  HEADTLAT  "NlTT  HEAD  (  LATS1  )  »  AHEAD  (M  )  , 
1   AMFAntLONr.CHK  ) 

W°TTF  OUTPUT  TAPF  JM ,  21,  (AHEAD(I),  I=M,  LONGCHK)  "~ 
fLl2LL-WR-LJA„0-UT?UT  TAPE  JL.28.  HEAD!  J)  .  J_BAR0Mm_,_I_C2  (  I  )  .  1  =  1.10). 

i "  t  ic5  m  »"i C3 m"»  i  ="i  »T6T»  fic4 1 1 ) .  ic3 1 1  > .  i  =1,1 0 1 " .  " 

__  2 I :  ( IWINDVFH I,L) ,  I C  6  (  I  )  ,  1  =  1.10).  L  =  J . J ) __ 

"401?  "tF  (MARKR)  ?25,226»225 

225  WRITE  OUTPUT  TAPE  JM.27,   HEAD(J),  (  XREJ  (  L  )  .  L=  1  .  20  ) 


?7f,  COMTINUF 
._      ™_„„„.„____„__„_„_.___ 

TABLE  (3  )_=J_  A R L_F _<_3  )+_ 0_._0 0_1  __ 

LINE"=0 
7?77  M  =  |.QNGCHIC/7 __ 

|_  =  M_4 
401 5  WRITE  OUTPUT  TAPF  JL»  19. 

WRITF  OUTPUT  TAPE  JL»  44,  TARLEm.  AMONTH.  IYEAR"" 

WRITE  OUTPUT  2AiiL.^i--!__?P-!.-Jl?A9lAi^J^.L_^_H_EAp2_(_L_ATS_2j  .AHEAD2  (  L)  , 

1    AHFAO?"(  M) 
WRITE  OUTPUT  TAPE  JL.  22.  ( AHEAD2 ( I ) .  I=L.M) 


4017  DO  229  J=LATN2.  LATS? 
J_P__5_OQ4___I=_l_,5 

"  Vi  1 1  )=jsFA<r.jrr~ 

IF  (ICKD)  228,226,227 

?  ?  7~  ir?  ( i )  =  jsf  ate'mp  ( i  ♦  j  i  / 1  c  if  1 1 

xir?(  n  =  ir?(n 


xir?t  n  =  xic?(  I  )/io.o 
_5Aj2QMiJLl=.LX.L£2JJ.J_-32.»0.iyjLJJl- 

RO  TO  5^04 
.__22J_-J_l__2__-__=-l-_- 

BAPOM( I )=0.0 
^0^4  rOMTTMMF 


LIME=LINF+1 

IF._U_Ji4EL-3_ril__5229_i.52.25.taQl* 

9014  TARLF(3)=TARLF(3)+O,001 
L-I_NE__1 

9<">15  WRITF  OUTPUT  TAPF  JL.  19. 
WPTTF  OUTPUT  TAPF  .11.  44.  TARLFI3).  AMONTH,  IYFAR 


WRITF  OUTPUT  TAPF  JL .  20,  HEAD2 ( L ATN2 ) ,  HEAD2 ( L ATS2 ) , AHEAD2 ( L ) , 
J ___EA__2jL__) 

WRITF  OUTPUT  TAPF  JL.  22.  (AHEAD2II).  I=L.M) 

-A2-?_9--WPaT___jaLiieiii__TAj-F__jj___^^  .. 


1  1  =  1. ■>) 
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_L2_---_--_UMJ__ 

r FTVF    nFGRFF    OUTPUT     LOOPS... 

TARLE(4)=TARLF(4)+0.001 

1_J__UL___. 

Un]  o    M_|_ONGrHK  /5 
L_=---l 

4"?i  V-ITF  OUTPUT  TAPF  JL  »  1 9 » 

I.IPITF  OUTPUT  TtPF   II.  65.  TARI  Fltl.  AMONTH.  IYFAR 

WRITE  OUTPUT  TAPE  JLi  51.  HEAD5 ( L AT N5 ) .  HE AD5 ( LATS5 ) . AHEAD5 ( L ) . 
1 AHEAD5-C.M.1. 

WR1TF  OUTPUT  TAPE  JL.  23.  (AHEAD5M).  I=L.  M) 


.AQZ2..DH.2.fiJh..A=LAlJi5j..LALS5.. 
HO  245  1=1,2 

TTFg,T  =  Q 


ITT ( I )=KSFA( I, J) 

J.E-jU£.1.LIJl1_23l3.j^.3.1«23£L 

23  0  XTG7  1  I )=KSFATFMP( I »J)/IC1 ( I  1 

JLLC2.1JL1=JSJ-C2.LI_LZJLQjlQ 

GO    TO    5^05 
73?     XTr?(I)=Q.^ 

!TFST=1 
5005    IC3tn=JDEWt  I,J1 

IF    (103(11)    235,235.234 
234    XIC4( I)=JDEWPT( I,J)/IC3(I ) 

x~i  c~4  <  i )  =Tx  i  C4  u )  -  3  2  6 . 6 1  / 1 8 .  n 

GO    TO    5006 


??<5    xTr4(I)=n,n 

_JJE__T_=_1 

Rn^A    ir5< I )=kwimD( I »j) 

I£_JLK51JJ-l-23.7.i23J-t2.3J'- 

7^6     TG6( I )=KWIMniRX( T .J) /IC> ( I ) 
Tr7in=KWTNOIRYl  I  .J)/IC5(  T) 


GO    TO    5007 
...  2-17„K.6JLI_L=il 

!07( I )-C 

ULESlal 

■5007     IC8(  t  )=KWINDV(  I  ,J) 
IF     (IC8(I>)     739,739,738 


?^o    XIC<»(  I  )=KWIN0VFL(  I  ,J)/IC8  (  I  ) 

GQ..XO-J anna. 

7^o   xiro(  n=o.n 
JJLEAI_=1 

snop    inot  I1=JCL0(  I  .J) 

TF  nriPIUI  741  .741  .740 


74  0  OLOUrX  I )=JCL0U0( I,J)/IC10(I) 

CLOUDi_J_=CiO_lDiI_U.L_-___L 

GO  TO  760 
-2A___Cl_0U0±_J-=_-_D_ 

ITFST=1 
76n  rr7m=KR6BiT..n 


IP  (ir?(!n  767. 762.761 
— J7-6.J— BJiBIiHJ.lXs.iUi.iS.QH.l-LtJ.U-lCZ.lXX 

pa  ROM! I )=RAR0M( T )/in.n 
Ga_iOL.5ja.QSL 

767  barom(  I  )  =  n.,ri 

ROOO  TF  irnilNTI  T..M  1  501  1  .5011  .5010 


5010  n|FF(I.J)=nlFF( I, J) /COUNT  (  I .J) 

__l_iO_-5_ll_- 

5011  DIFFI I .J)=0.0 

IIESIaI 

501?  IF  (XINTF(COUNT( I .J) 1+IC1 ( I ) )  290.290,291 
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791     A T otfmp ( n  =X ! C2( ! ) -0 [FF( I , J ) 

fin  Tfi~?92~ 

?if    A  tptFmo(  I  )  -o.o 

79  7    KKi<l<:  =  KKk:K  +  x^AXOF(IC1<I),IC3<I>»!C5<I>,IC8(I>.ICl(JlI>.IC2(I>) 
IF    (ITEST)    5014,5014,5013 

"Vni  i   o>"Tl  f=o  . n 

-  01 1 "  =0_tp 

OF(I)=o.n 

0"tT)=n.o 

0H( I )=0.0 
GO    TO    ?45 

5014    if    <XMINOF(  IO  (  I  )  .IC3<  I  )  .IC8(  I  )  .IC101  I  H-5)    5ol3»5o~15  t5ol5 
.^J^_XEXJr?j.lJ_+J_to. 

I  =fXTC7(T  )-X )»1".0  +  1.p 

M=XTC4(  I  )+l .0 
JC=M=1 

N=(XIC4( T )-X)»1o.o+i .0 
X=FS(L,K)-FS(N,M) 

QR( I )=(  1.14E-07)*(XIC2  (  I  r+273.dl**A*l0.39-U.U50*EMN»M)"**U.5)* 
1       ( 1 .  0-XK2( J)»CLOUD( I)»»2)+<4«56E-07)*(XIC2t I ) +273.0 ) ««3«DI FF( I .J) 


qo(  I  )=-0R( I ) 
QT( T  )=HFAT(N0,J)»(1.0-X1C1 ( J)*CLOUD( I ) -Q.38*CL0uDt 1 )»»2) 

0R(  I)=QI(  I  )*RFFLFCT(NO»J)*(~-1.0) 
0F(Tl=?.»1»XTr9(T  )*X»t-1  .0) 

7  4  3    0H(  I  )=-1  .*8*r>IFF(  I  ,  J)*XICQ<  I  ) 

76  5    nTiriifiiii  )+op>  r  n-nst  t  uo^i  i  i+qh(  t  i 


LINF=LINF+1 
LF..lJ_I^.Er.6J..52&Aji2iL6jL9ni7- 

0017  TAPLE(4)=TAALF(4)+0.001 

J_I.!AE=1 

9(118  WRITF  OUTPUT  TAPF  JL »  19. 
WRTTF  OUTPUT  TAPF  JL,  45.  TARLE(4).  AMONTH.  IYEAR 


M=L0NGCHK/5 
J_=M=d. 

WRITF  OUTPUT  TAPF  JL.  51.  HEAD5 ( LATN5 ) .  HE AD5 ( L A  I S5 ) . AHE AD5 ( L ) . 

J AU£AI151^J 

WOITF  OUTPUT  TAPF  JL .  23.  ( AHEAD5 ( I ) »  I=L»  M) 
S746  WRTTF  OUTPUT  TAPF  JL.30.  HEAP5(  J  )  .  i  X  IC2  1  L  )  ,  I  CH  L  )  ,L=  1  .  2  )  > 


1    (AIRTFMP(L).  TCKL),  L  =  1.2).  ( X IC4( L ) . I C3 ( L ) > L= 1 »2 ) » 
?(BAROM(L) ,TC7(L),L=1 ,2)  ,  ( IC6( L) . IC5 (L) ,L=1>2) » I IC7I L ) » IC5 1 L ) » 

3    L=1.2).(XIC9(L).IC8(L),L=1.2)» (  CLOUD  <-L) »IC101L).L=1»2) 
J^RJ_LE_^yj_Py_T__TA_P_E_.JL.4.8_. (Ql 1L) ,QI  (  L  )  »QR  \  L)  ,QB(  L  )  .QElL)  > 

]    0H(L)»L=1»2) 
746  rOMTINUF 


c   ♦  *»*»♦*#»*»#«♦**********#**********«*** 
5747  IF  1  ALPHA ( 71 )-5H**»»*)  6969,248,6969 

r      PULL  PRFCFDING  CARP  RETURM  FOLLOWING  CARD 

(-5967  TF  (Al  PHA(71  )-5H»»»*»)  747.748,247 

747  LOMfirHK,=LOMGrH*'  +  10 
TF  M  ONGCHK-I  IMI  ONGE )1 035 . 1 03 5 ♦ 248 


C BUFFER  ERROR  ROUTINE 

7000  BACKSPACE  JK 
ICnUJiI=_lC_QUKI_-tl 

IF  (ICOUNT-5)  7001  .7001 .7002 
70Q1  RUFFFR  TN  l.X.Ol  (  T PUFFFR I<  1  )  . I RUEFFR 1 ( 1 5 )  ) 


»#*»*»♦#» 


GO    TO    NOT,     (6000,6002) 
.JLQJl2-J3Ll£F£B._liJ„J^iKj.0X-_LiaUEEFji1_(.i.uJLBU£E£Bj-il5JLl. 

GO    TO    7nnl 
-aOn3__5F4i(4-MD— JJC 

jtc  =  jk:  +  i 
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.__GQ„rO__7£QJ _       . 


?4*     LLLL  =  ^^  K'K-t-JJJJ 


WRITF    OUTPUT    TAPE    JL  .    32.     AMONTH , I  YEAR , XDATE , I  DAY ♦ JYEAR 
280    W_RITF    OUTPUT    TAPE    JL.47.LI_LL  i  JJJJ»KKKK»XNA_ME    _ 
"7^r^l?ITE"~6UTTuT~~t~AP¥"LT4¥»^M6NTHt  1  YEAR  .TABLE  1  1  T.AMONl  H.  I  YEAR. 

1        T  ARLF  (  2  )  .  AMONTH.  I  YEAR .TABLE ( 3 ) 

~"wRITe~0UtP~UT    Ta'Pe"  Jl",50.  AMONTH.  I  YEAR.  I  ABLeTV) 
771     FMD    F!LF    Jl 


JF     (MARKT+MARKR)     249.250.749 

7frP       FMP|      F  T1    F      .»M 

75n     e;TOD 

FMn 

TOFMT            ~~ZFPO   " 
FMfPY ZF^O 


7rrjn  pyj  *# 

iAii l±. 

LLS  48 
_S.AU. Z2_ 

Z.  LIL  «j      ** 

FMA  <"> 


7"»  TJP  ft       *  +  ! 

AUL ZJL51L 

7?                        STA              «,        ** 
J_L_J Z3__ 

FND 
TDENT HOLT 


FMT"Y      DOLT 

_V___t __X__ ** 

SAU        DOLT8 

LLS. 48 

SAU        *+l 

LIL     ft 


»* 


IN!      ft   -1 

_50|_Tg __.___>. ** 

FMA        n 

LQ0-°_ 1±5l 6 

TNA        -12R 
AJP     7       *  +  ? 


AJP     P   DOLT 

_fm_a_ n 

! JP     6   LOOP 
__  _J_LJ DOLT. 

FND 


FMD  FILE  T 
_Q.FW_J_Nn.__J 

T  =  T  +  1 
_X_=X_+1_5___T_ 

P FT URN 

FMD 


FMD        PUCOM 
FINIS 


-FXFTUTFR. 
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APPENDIX  B 


Tables  used  in  data  processing  program. 
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Appendix  Table  B-l. — Values  by  month  and  by  latitude  of  incoming  solar 
radiation  with  a  clear  sky.   [Adapted  from  Berliand,  I960. ] 


JAN 


FEB 


MAR 


APR 


MAY 


JUN 


JUL 


AUG 


SEP 


OCT 


NOV 


DEC 


55N 

80 

173 

355 

548 

696 

761 

713 

570 

398 

229 

109 

58 

50N 

131 

237 

411 

589 

718 

774 

732 

609 

450 

288 

162 

105 

45N 

190 

305 

464 

625 

738 

784 

752 

645 

497 

348 

223 

162 

40N 

255 

371 

514 

655 

753 

789 

767 

675 

539 

405 

287 

227 

35N 

321 

431 

559 

679 

762 

789 

774 

697 

578 

458 

347 

290 

30N 

381 

485 

597 

696 

764 

784 

773 

711 

612 

507 

403 

349 

25N 

437 

531 

62  7 

707 

759 

774 

766 

718 

641 

551 

456 

405 

20N 

487 

571 

652 

710 

747 

759 

752 

718 

663 

590 

506 

458 

15N 

533 

607 

672 

708 

731 

737 

732 

713 

681 

622 

551 

507 

ION 

575 

637 

688 

702 

709 

708 

708 

704 

693 

649 

591 

553 

5N 

615 

661 

700 

693 

682 

674 

678 

689 

700 

671 

627 

597 

ON 

651 

680 

706 

680 

651 

637 

645 

670 

700 

689 

658 

637 

OS 

681 

696 

706 

663 

619 

599 

609 

647 

693 

700 

685 

675 

5S 

709 

710 

699 

643 

585 

558 

570 

618 

679 

705 

708 

710 

10s 

733 

719 

687 

616 

547 

513 

529 

585 

659 

703 

726 

740 

15S 

753 

724 

670 

584 

506 

467 

486 

550 

635 

695 

740 

764 

20S 

769 

726 

650 

547 

462 

419 

440 

511 

571 

685 

751 

784 

25S 

782 

722 

625 

508 

416 

372 

391 

469 

539 

671 

758 

802 

30S 

790 

710 

595 

468 

368 

323 

341 

423 

540 

652 

760 

816 

35S 

792 

691 

559 

424 

315 

269 

291 

374 

502 

629 

756 

826 

40S 

791 

667 

520 

377 

260 

211 

236 

324 

461 

600 

74  5 

831 

45S 

784 

638 

477 

328 

206 

153 

178 

269 

416 

566 

729 

828 

50S 

771 

604 

430 

274 

152 

100 

124 

211 

367 

527 

711 

818 

55S 

753 

567 

380 

215 

102 

54 

75 

152 

310 

486 

696 

808 
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Appendix  Table   B-2. — Values   by  month  and   by   latitude  of   the   percent  of    incoming 
solar  radiation  that    is   reflected.      [Adapted   from  Budyko,    1956.] 


JAN  FEB  MAR  APR  MAY  JUN  JUL  AUG  SEP  OCT  NOV  DEC 

55N  .190  .150  .105  .077  .077  .067  .077  .085  .095  .132  .177  .197 

50N  .170  .130  .095  .072  .072  .062  .072  .075  .085  .117  .152  .172 

45N  .147  .112  .087  .067  .067  ,060  .067  .067  .077  .102  .132  .150 

40N  .122  .097  .082  .062  .062  .060  .062  .062  .072  .087  .117  .130 

35N  .105  .087  .077  .060  .060  .060  .060  .060  .067  .077  .102  .112 

30N  .095  .082  .072  .060  .060  .060  .060  .060  .062  .072  .087  .097 

25N  *085  .077  .067  .060  .060  .060  .060  .060  .060  .067  .077  .085 

20N  »075  .072  .062  .060  .060  .060  .060  .060  .060  .062  .072  .075 

15N  «067  .067  .060  .060  .060  .060  .060  .060  .060  .060  .067  .070 

ION  «062  .062  .060  .060  .060  .060  .060  .060  .060  .060  .062  .070 

5N  .060  .060  .060  .060  .060  .060  .060  .060  .060  .060  .060  .067 

ON  «060  .060  .060  .060  .060  .060  .060  .060  .060  .060  .060  .062 

OS  .060  .060  .060  .060  .060  .062  .060  .060  .060  .060  .060  .060 

5S  .060  .060  .060  .060  .060  .067  .060  .060  .060  .060  .060  .060 

10S  .060  .060  .060  .060  .062  .070  .062  .062  .060  .060  .060  .060 

15S  .060  .060  .060  .060  .067  .070  .067  ,067  .060  .060  .060  .060 

20S  .060  .060  .060  .062  .072  .075  .075  .072  .062  .060  .060  .060 

25S  .060  .060  .060  .067  .077  .085  .085  .077  .067  .060  .060  .060 

30S  .060  .060  .062  .072  .087  .097  ,095  .082  .072  .060  .060  .060 

35S  .060  .060  .067  .077  .102  .112  ,105  .087  .077  .060  .060  .060 

40S  .062  .062  .072  .087  .117  .130  .122  .097  .082  .062  .062  .060 

45S  .067  .067  .077  .102  .132  .150  .1*7  ,112  ,087  .067  .067  .060 

50S  .072  .075  .085  .117  .152  .172  .170  .130  .095  .072  .072  .062 

55S  .077  .085  .095  .132  .177  .197  ,190  ,150  ,105  .077  .077  .067 
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Appendix  Table  B-3. --Values  for  the  constant  k  by  latitude  used 
in  computation  of  back  radiation.   [Adapted  from  Budyko,  1956.] 


55N  50N  45N  40M  35N  30N  25N  20N 

0.75  0.73  0.71  0.69  0.67  0.64  0.62  0.60 

15N  ION  5N  ON  OS  5S  10S  15S 

0.58  0.56  0.53  0.51  0.51  0.53  0.56  0.58 

20S  25s  30s  35S  40S  45s  50s  55S 

0.60  0.62  0.64  0.67  0.69  0.71  0.73  0.75 
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Appendix  Table  b-k. —  Values  of  saturation  vapor  pressure  over  water  used  in  the 
computation  of  evaporation.   [Adapted  from  Smithsonian  Meteorological 
Tables,  1958,  Table  3k,    p.  352.] 
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APPENDIX  C 

Sample   page  copies  of  data  tabulations  produced  by  the 
FORTRAN-62  data  processing  program  listed  in  Appendix  A. 


Note:    Computer  printout  labels  and  sign  convention  in  Table  4.022  differ 
from  the  text  (see  p.  17)  as  follows: 

1) 


Q(T) 

=   QT. 

Q(D 

=  Qr 

Q(R) 

=  %> 

Q(B) 

-  Qb- 

Q(E) 

=  QE- 

and  Q(H) 

-  %■ 

2)      Q(T)  =  Q(I)  +  [-Q(R)]  +  [-Q(B)]  +  [-Q(E)]  +  [-Q(H)] 


This  equation  is  equivalent  to  that  given  on  page  17.  The  signs  in  the 
equation  and  computer  output  values  have  been  altered  to  permit  proper 
notation  of  the  components  (heat  loss  from  sea — negative;  heat  gain  by 
sea — positive). 
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U.S.    BUREAU  OF  COMMERCIAL  FISHERIES 

BIOLOGICAL  LABORATORY 

SAN   DIEGO,  CALIF. 

MARINE  WEATHER  OBSERVATION  SUMMARY 

FOR  THE 

PACIFIC  OCEAN 

FOR    THE  MONTH  OF 
JUNE    1964 


DATE  OF  RUN JUNE  30,   1964 


NUMBER  OF  OBSERVATIONS  9241 

NUMBER  OF  REJECTIONS  391 

SUMMARIZED  OBSERVATIONS        8850 


PROGRAMMED  BY.  .  .      MARVIN  W.  CLINE 
PROCESSED  BY JAMES  A.  RENNER 
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Created  in  1849  the  Department  of  the  Interior — a  depart- 
ment of  conservation — is  concerned  with  the  management, 
conservation,  and  development  of  the  Nation's  water,  fish, 
wildlife,  mineral,  forest,  and  park  and  recreational  re- 
sources. It  also  has  major  responsibilities  for  Indian  and 
Territorial  affairs. 

As  the  Nation's  principal  conservation  agency,  the  De- 
partment works  to  assure  that  nonrenewable  resources  are 
developed  and  used  wisely,  that  park  and  recreational  re- 
sources are  conserved  for  the  future,  and  that  renewable 
resources  make  their  full  contribution  to  the  progress,  pros- 
perity, and  security  of  the  United  States — now  and  in  the 
future. 
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